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ABSTRACT 

Collision  integrals  developed  and  partially  evaluated  by 
L.  W.  Hunter  and  R.  F.  Snider,  and  utilised  by  Hunter  in  his 
treatment  of  the  effects  of  a  magnetic  field  on  the  shear 
viscosity  and  thermal  conductivity  cf  single  component  diatomic 
gases,  are  cast  in  the  Generated  Phase  Shift  (GPS)  formalism 
of  C.  F.  Curtiss.  This,  along  with  the  introduction  of  certain 
operators,  allows  the  collision  integrals  to  be  considerably 
simplified  by  facilitating  the  evaluation  of  several  summations 
and  angle  integrations. 

The  difficulties  inherent  in  treating  diatom-diatom 
interactions  lead  to  consideration  of  binary  gaseous  mixtures 
in  which  the  dominant  species  is  atomic,  the  diatomic  species 
being  restricted  to  low  concentrations,  'inch  systems  require 
consideration  of  atom-atom  and  a  tom  di.tf'T  interactions  only. 
Employing  methods  introduced  by  Hunter,  scalar  equations  are 
obtained  for  the  transport  properties  of  binary  mixtures  in 
applied  magnetic  fields.  The  collision  integrals  occurring  in 

♦This  research  was  supported  by  National  Vim,  e  foundation 
Grants  CHE 74-1  /4'J4  AO!  and  CHE77- I/i/f.l  . 


these  equations  are  found  to  be  generalizations  of  those 
discussed  by  Hunter  and  Snider. 

The  shear  viscosity  of  an  atom-diatom  mixture  in  an  applied 
magnetic  field  is  then  treated  in  detail.  The  basis  set  is 
truncated  and  the  diatomic  species  is  restricted  to  low 
concentrations.  The  expressions  obtained  for  the  shear 
viscosity  tensor  are  in  qualitative  agreement  with  experimental 
observations.  Steps  leading  to  calculation  of  the  shear 
viscosity  tensor  in  this  particular  case  are  discussed. 
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IN'L’ODUaiON 


1.  The  Senftl eben-Beenakker  effect 

Calculated  values  of  the  transport  properties  of  dilute 
polyatomic  gases  that  compare  favorably  with  experimentally 
measured  values  can  be  obtained  by  treating  the  molecules 
as  point  particles,  i.e.,  as  structureless  particles 
interacting  via  a  spherical  intermolecular  potential 
Incorporating  internal  molecular  structure  and  angular 
dependence  of  the  interaction  potential  complicates  the 
transport  property  calculations  considerably,  but  makes 
only  small  contributions  to  the  calculated  value.  This 
fact  somewhat  reduces  the  incentive  for  carrying  out  such 
calculations,  and  makes  it  difficult  to  assess  the  relative 
merits  of  various  treatments  of  the  non-soherical  contribution. 
The  spherical  treatment,  however,  is  incapable  of  accounting 
for  the  effect  of  an  anplied  external  field  on  the  transport 
properties  of  polyatomic  molecules.  This  effect,  known  as 
the  Senftl eben-Beenakker  effect  ♦ jr  too  men  who  first  observed 
it  experimental  ly  for  pararaure  f  i.~  and  diamagnetic  molecules, 
respectively,  can  actually  be  measured  more  accurately  than 

O 

the  transport  properties  ttiemsol  ves  .  It  is  extremely 
sensitive  to  tiie  anisotropic  part  of  the  intermolecular  potential 
of  molecules  possessing  internal  structure,  and  for  this  reason, 
may  provide  a  useful  means  of  probing  nonspherical  interactions. 


There  is  a  simple  physical  explanation  of  the  effect. 
Under  the  influence  of  a  gradient  in  one  or  more  of  the 
macroscopic  variables  of  a  system,  the  molecules  tend  to  change 
their  distribution  so  as  to  reduce  tne  magnitude  of  the  gradient, 
or  gradients  involved.  The  result  is  a  flow  within  the  system. 
Since  the  molecules  are  nonspherical ,  the  collision  probabilities 
are  functions  of  the  orientations,  and  consequently,  the 
molecules  have  a  tendency  to  align.  Inis  results  in  a  value 
for  the  transport  property  that  is  slightly  larger  than  that 
which  would  obtain  in  the  absence  of  any  alignment. 

Since  both  paramagnetic  and  diamagnetic  molecules  have 

net  maqnetic  moments,  in  the  presence  of  a  magnetic  field 

they  precess  about  the  field  direction,  except  for  those 

components  parallel  to  the  field,  this  precession  tends  to 

destroy  the  preferential  alignments,  or  polarizations,  resulting 

from  the  presence  of  the  gradients.  As  a  consequence,  the 

transport  properties  of  molecules  with  net  magnetic,  moments 
a 

are  reduced  in  the  presence  of  an  externa)  magnetic  field. 

The  extent  of  the  destruction  of  tne  polarizations  is 
directly  related  to  the  average  number  of  precessions  between 
collisions.  The  precession  frequency  is  proportional  to  the 
magnetic  field  strength,  ,  while  the  time  between  collisions 
is  inversely  proportional  to  the  pressure,  p  ,  for  a  dilute 
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gas.  At  constant  temperature,  the  extent  of  the  destruction 
of  polarization,  and  consequently  the  Senf tleten-beenakker  effect, 
is  solely  a  function  of  Hf,/p  •  At  sufficiently  high  field 
strengths,  polarization  is  cor  pletely  destroyed,  and  saturation 
occurs . 

In  extending  the  above  picture,  it  is  apparent  that  any 
polyatomic  gas  will  exhibit  a  similar  effect  when  subjected  to 
an  external  field  capable  of  causing  the  molecules  to  precess. 
Other  than  the  case  of  polar  molecules  in  electric  fields, 
however,  the  field  strengths  required  to  cause  a  reasonable 
degree  of  precession  are  so  high  as  to  make  practical 
applications  virtually  impossible. 


2.  Theory  and  Experiment  -  A  Grief  History 


t  ,  ’) 


The  effect  of  a  magnetic  Mold  the  thermo !  conductivity 

of  gaseous  02  was  first  observed  by  h.  Sonftleber,  ‘  in  1 030 . 

Engelhardt  and  Sack ^  then  observed  the  Senftlenen  effect  on 

the  viscosity  of  02  in  1932.  ’he  f  act  teat  NO  bor.a  v s  1  r  the 

same  manner  as  02  coupled  with  trie  fact  that  in  matures  with 

non-paramagnetic  gases,  the  efv<  t  is  r^noor tion a  i  to  the  role 

fraction  of  the  paramagnetic  species,  led  early  observers  to 

the  (erroneous)  conclusion  that  the  Senftleben  effect  is  a 

property  of  paramagnetic  gases  only.  Other  observations  o* 

a 

the  first  decade  of  sucn  studies  were  that  in  trie  presence 
of  a  magnetic  field,  the  transport  coefficients  of  paramagnetic 
gases  decrease  by  from  0.1  to  1.0':;  that  the  effect  is  even 
in  the  magnetic  field;  and  that  at  constant  temperature,  it  is 
a  function  of  the  field  strength  divided  by  the  pressure. 

The  observations  of  the  first  decade  led  C.  J.  Gorter ^ 
to  an  early  qualitative  explanation  of  the  effect  in  1938. 

His  explanation,  based  upon  con..  ide*~at  i  :n  r.f  tec  effect  of  a 
magnetic  field  on  the  mean  free  path  o*  an  j,  molecule,  was 

o 

treated  more  quantitatively  in  1939  by  c’ernUe  and  van  Lier. 


These  early  developments  treated  the  rotating  paramagnetic 
molecule  as  a  disc  with  a  magnetic  moment  in  the  direction  of 
the  axis  of  rotation.  The  cross  section  depends  on  the  angle 
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between  the  axis  of  rotation  and  the  direction  of  motion.  In 
the  absence  of  a  magnetic  field,  the  direction  of  the  axis  of 
rotation  is  conserved  between  collision:,.  In  the  presence  of 
a  field,  however,  the  magnetic  moment,  and  consequently  the 
axis  of  rotation,  precesses  about  the  direction  of  the  field. 

The  collision  cross  section  now  changes  periodically  during 
the  flight  of  the  molecule,  necessitating  an  averaging  over 
the  precessions  in  the  mean-free-path  picture.  Tin's  additional 
averaging  leads  to  the  decrease  in  the  transport  properties. 

The  mean-free-path  approach  adequately  accounts  for  the 
phenomenon  of  saturation,  the  functional  dependence  on  HM/p  , 
and  the  dependence  on  the  mole  fraction  of  pu ; amagnetic  species 
in  a  mixture. 

Following  the  1930's,  during  which  the  Senftleben  effect  on 

paramagnetic  gases  was  studied  extensively,  essentially  no 

studies  of  field  effects  took  place  until  tne  early  1960's. 

i  o 

Then  in  1961,  Kagan  and  Maksimov  '  r^f-senl'mi  a  variational 
calculation  of  tne  Senftleoen  effect  **w  v  o't.jI  conductivity 
of  Or.  Their  treatment  of  the  1  ”  •  •-  /-d  'cm  ca  1  Pcltzmann 
equation  utilizes  an  elastic  coll  is  on  model  and  a  truncated 
expansion  of  the  di  stri  buticn  funct’on  in  terms  of  irreducible 
Cartesian  tensors  in  the  (reduced  peculiar)  velocity,  W  ,  and 
the  rotational  angular  momentum,  ,j  . 


As  early  as  1939,  Chapman  and  Cowling  had  realized  that 
in  the  presence  of  a  perturbing  gradient,  for  exjmpl**  1  t.t  T  , 
the  nonequilibrium  distribution  unction,  wri » ten  as 
f  =  f^(l  +  t> )  with  .*  =  -A-V  i,:  T  ,  could  contain  terms  involving 
J  .  Since  no  physical  reason  .vas  seen  to  consider  polarization 
of  nonspnerical  molecules  in  dilute  gases  at  tnat  tin*.*,  they 
disregarded  the  anisotropy  of  the  distribution  function  in  .1  , 

choosing  instead  to  write  A  in  the  monatomic  form,  A(W-')W  . 

1 2 

In  1961,  in  a  treatment  of  diatomics,  Kagan  and  Afanas’ev 
noted  correctly  that  A  should  include  terms  involving  all 
possible  true  vectors  which  can  hr*  constructed  from  W  and  .!  , 
i.e.,that  A  should  be  written  as  A  =  A  t  VI  *-  A  :W  *  J  *"A-j(W*_J).J  . 

The  Senftleben  effect  is  then  the  result  of  the  reduction  of  the 
anisotropy  in  j  space  due  to  the  precession  of  the  molecules 
in  the  presence  of  a  magnetic  field. 

13 

The  following  year,  1962,  J.  J.  M.  Beenakker  measured  the 
change  in  the  viscosity  of  N9  and  other  diamagnetic  gases  in  the 
presence  of  a  magnetic  field,  removing  the  restriction  of  these 
field  effects  to  paramagnetic  molecules.  Analogous  measurements 
of  the  thermal  conductivity  were  conducted  by  Gorelik  and 
Sinitsyn  ^  shortly  thereafter.  This  experimental  work  involving 
field  effects  on  diamagnetic  gases  has  been  followed  by 
considerably  renewed  interest  in  the  Senf tleben-Beenakker  effect 

5 

from  both  experimental  and  theoretical  standpoints.  Extensive 


experimental  information  is  now  available,  as  well  as  numerous 
theoretical  developments.  Unfortunately,  very  few  numerical 
calculations  of  field  effects  on  transport  properties  have  been 
attempted.  This  work  is  a  step  toward  such  calculations. 
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3.  Scope 

To  date,  most  theoretical  development.1  af  the  rransport 
properties  of  polyatomic  qa.es  in  applied  fields  utilize  a 
Chapman-Enskog  procedure  ^  to  solve  the  equations  for  the 
transport  coefficients.  Arguments  are  advanced  to  justify  the 

truncations  made,  but  most  are  not  truly  justified  until  the 

2 

solution  has  been  ob;.v-.ed.  L.  W.  Hunter  hr,  devised  a 
procedure  to  treat  transport  properties  that  alio, vs  the  tensor 
analysis  to  be  completed  before  any  truncation  is  i.o  rssary. 
Hunter's  techniques  are  examined  in  "art  f  of  th i  .  thesis ,  and 

the  collision  integrals  arising  in  his  treatment,  ao-  cast  in 

I  S 

terms  of  the  reduced  scattering  matrix  or  Curtiss  and 
co-workers.  These  collision  integrals  are  then  simplified 
considerably  by  the  introduction  of  certain  operators  and  the 
explicit  evaluation  of  a  number  of  the  sums  and  angle 
integrals  involved. 

Hunter's  treatment  applies  to  a  s ingle-component  diatomic 
gas.  Calculation  of  the  cr, import  nro.mvt  ies  or  such  a  system 
is  complicated  by  the  dependorv  e  of  the  interaction  potential  on 
the  orientations  of  the  colliding  diatomirs.  Recently,  a 
calculation  of  the  transport  properties  of  an  atom-diatom  mixture, 
in  which  the  diatom  is  present  in  low  concentration,  has  been 
completed  by  R.  Wood.  ^  Such  a  calculation  requires 
consideration  of  atom-atom  and  atom-diatom  interactions 


exclusively.  Since  the  Senftleben-Beenakker  effect  is  observed  in 
mixtures,  it  seems  reasonable  that  a  detailed  calculation  of  the 
effect  should  first  be  attempted  for  an  atom-diatom  mixture 
involving  a  diamagnetic  diatomic  species  in  low  concentration,  and 
later  generalized  to  a  purely  diatomic  gas.  In  establishing  the 
groundwork  essential  to  an  atom-diatom  calculation.  Hunter's 
development  for  a  single  component  gas  is  first  generalized  to 
binary  gas  mixtures.  This  is  accomplished  in  Part  II. 

Then  in  Part  III,  the  equations  derived  in  Part  II  for  a 
binary  gas  mixture  are  restricted  to  an  atom-diatom  mixture  that 
is  predominantly  atomic,  and  expressions  for  the  shear  viscosity 
are  developed  in  detail.  The  expressions  obtained  are  shown  to 
be  consistent  with  experimental  results,  and  the  collision 
integrals  involved  are  expressed  in  such  fashion  as  to  make 
evident  the  extension  of  the  calculations  of  R.  Wood  to  the 


present  problem. 
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A  SINGLE  COMPONENT  DIATOMIC  GAS 
IN  AN  APPLIED  MAGNETIC  FIELD 


1.1  Transport  Coefficients 

2  3 

The  linearized  Wal dmann-Snider  equation  ’  can  be  used  to 
obtain  a  set  of  tensor  equations  from  which  the  transport 
coefficients  of  a  dilute  gas  can  be  determined.  In  most 
appl ications ,  the  number  of  polarizations  of  the  gas  is  limited 
to  a  finite  set  and  the  tensor  analysis  performed  in  a  truncated 
basis.  L.  W.  Hunter  has  developed  a  technique^  for  carrying  out 
the  tensor  analysis  with  the  complete  set  of  polarizations,  i.e. 
before  truncation  of  trie  basis.  In  Reference  1,  Hunter  demonstrates 
this  technique  for  the  linear  Zeeman  effects  on  the  shear  viscosity 
and  thermal  conductivity  of  diamagnetic  molecules.  A  sketch  of 
Hunter's  derivation  of  the  scalar  transport,  equations  is  given  in 
this  section.  The  analogous  development  for  binary  mixtures  is 
carried  out  in  detail  in  Part  II. 

In  Hunter's  work,  the  equations  for  the  transport  coefficients 

involve  collision  integrals  which  are  related  to  those  discussed 

4 

and  partially  evaluated  by  Chen,  Moraal ,  and  Snider  and  by 
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5 

Hunter  and  Snider.  For  convenience,  expressions  hoi;;  Reference  1 

will  be  referred  to  by  their  equation  number  (from  the  reference) 

preceded  by  H,  e.g.  H-40.  Similarly,  those  from  Reference  4  will 

be  preceded  by  CMS,  and  those  from  Reference  5  by  US. 

To  obtain  the  basic  tensor  equations  for  the  transport 

coefficients,  the  Wigner  distribution  function-density  operator, h 

/  —  \ 

f  ,  is  expanded  ahout  the  local  equilibrium,  f' 

f  =  f^0)(i  -  rf  + ...).  *i.i-i :> 

The  perturbation,  d  ,  is  written  in  terms  of  quantities  linear 
in  gradients  of  the  macroscopic  variables.  If  only  shear  viscosity 
and  thermal  conductivity  are  of  interest,1 

4>  =  -A*V  In  T  -  B:[7y0]^  ,  (1.1-2) 

in  which  the  bracket  notation  indicates  a  symmetric  traceless 
tensor,  T  is  the  temperature,  and  v0  the  stream  velocity. 

When  (1.1-2)  is  used  in  the  formal  expressions  for  the  heat 
flux,  q,  and  shear  pressure  tensor,  n ,  and  comparison  is  made  with 
the  phenomenological  expressions  for  these  quantities,  H-5  and 
H-9  are  obtained  for  the  thermal  conductivity  and  shear  viscosity 
tensors,  respectively: 

2  k  V 

*  =  (-jj^)  2n«(kBTW2  *  Hint)W|A»  (1.1-3) 


1 


D  =  nkRT  <-<[W]u;|B»  . 


(1.1-4) 


In  the  above,  kg  is  the  Boltzmann  constant,  m  is  the  mass  of  a 
gas  molecule,  W  =  (2mkgT)  2  p  is  the  reduced  peculiar  velocity, 
H.  is  the  internal  Hamiltonian  of  a  single  molecule,  and  the 
inner  product  is  defined  as  in  H-6: 


|b»  =  n’1  Tr  j  dp  f^  g  +  b  ,  (1.1-5) 


in  which  denotes  the  adjoint  operator,  and  the  trace  is  over 
internal  states. 

Due  to  the  independence  of  macroscopic  gradients,  the 
linearized  Wal dmann-Snider  equation  becomes  separate  equations^ 
for  A  and  B  : 


(r>.(  +  i  1 )  A 


1  m  | 


,  H.  .  -  <H.  .>  \ 

I  *  “ 


(1.1-6) 


(«  +  i L )B  = 


(1-1-7) 


in  which  or  t>  =  <J-'H.nt|l»  ,  d?  is  the  linearized  Waldmann-Snider 
collision  superoperator,  and  L  is  the  Larmor  precession 
superoperator  defined  by  H- 13, 


% 
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L  3  C  d  ^  *  r  3 


(1.1-8) 


In  the  above,  is  the  component  of  internal  angular  momentum 

operator  along  the  field,  and  is  the  Larmor  frequency: 


=  h  HM  ’ 


1.1-5) 


where  g  is  the  rotational  Landd  g-factor,  q  the  nuclear 
magneton,  and  the  magnetic  field  strength. 

Using  the  tensor  equations  just  obtained,  aion*  w i in  the 
auxiliary  condition  of  H-14, 


<CW|A»  =  0  , 


(1.1-10) 


scalar  equations  are  obtained  for  A  and  n  by  considering 

4  5 

two  bases,  the  projection  operator  basis,  ’ 


A 

=pqsj 


Lps(w)  v(q)(j)  P.(Pj)'y=, 


(i. i-i 


in  which  Lps(W)  is  a  p-th  rank  tensor- va 1 ued  function  of  W  , 

V^(J)  is  a  q-th  rank  tensor  operator  dependent  on  the 

internal  angular  momentum  operator,  J  ,  of  one  molecule,  P. 

is  a  projection  operator  onto  the  j-th  energy  level,  and  p^  is 

the  Boltzmann  weight  of  the  j-th  level;  and  the  Wang 

4 

Chang-Uhlenbeck  basis. 
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Bpqst  =  LPS(H)tJ](q)  ^’(^l  •  d.1-12) 

in  which  is  a  normalized  Wang  Chang-Uhlenbeck  polynomial 

and  [J]^  is  a  q-th  rank  symmetric  traceless  tensor  of  J  ,  e.g. 
[J]  =  ^  -  V3J  :'<J  ,  where  y  is  the  second  rank  unit  tensor. 

The  matrix  elements  of  R  and  L  in  these  two  bases  are  given 
in  H -24,  H-25 ,  H-40,  and  H-4 1 : 


<cA(a)a|fi:|A(f'Ja; .,» 

pqsj 1  1  pqsj 


(1.1-13) 


=  n(8kjT>’?''  (2at,r'  -(pqsj  «<<MV  . 


pqsj  1  pqsj 


Lpq  ^aa  *  6pqs j , p ’ q * s ‘ j ’  6u,a‘  ’ 


(1.1-14) 


«B(d)c>;|B(a'  j  , 
pqst1  1  pqst 


"'sSt)'1-'  (Za+l>'' 


pqst 
ip'q's't'j 


(a) 


6aa,a'al ’ 
(1.1-15) 


and 


pqst1  1  pqst' 


I  tA^(aa‘  )  d  t  <S 

pq  pqst,p  q  St  u,a 


(1.1-16) 


In  the  above,  .  is  the  reduced  mass  of  the  collision  pair,  and 

and  are  spnerical  components  of  the  total 

pqsj  pq  t 
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polarization  bases  of  A  and  6  .  ,  respectively. 

K  ,pqsj  -pqst 

Expressions  for  t^^(aa')  are  given  in  Ref.  1. 

The  spherical  components  of  A  and  B  are  then  expanded 

/  i  \ 

over  the  total  polarization  basis,  to  obtain 


A 


m 


y  ( i  -  <s 

pqsta 


pqs t , 1 000 


)  A1 


pqst 
( a  )m 


B(a)in 

pqst 


(1.1-17) 


(the  1-iS  s^.  ^qqq  factor  ensuring  that  the  auxiliary  condition 
expressed  in  equation  ( 1 . 1  - 1 0 )  is  satisfied)  and 


V  Rpqst  R(a)m 

pjta  f“'st  ’ 


(1.1-18) 


where  A?^1  and  are  expansion  coefficients  given  by 

( a )  m  ( a )  m 

the  following  matrix  equations: 


n(^T)'V?  (2a+l) 


-1 


fp  q  s  t 


I  ° 

p’q’s't'  p'q's’t' 


tP'q's'f 

(a)m 


+  i  ( 1  -  6 


pqst, 1000  L,  pq 


l  L(m)(aa’)APq^ 

L  nn  '  '  (a  )m 


fCintTl 


5kB  I  '! 


'■L  r  r  r  _  ' _ D__l  r 

3m  j  a,l L  pqs t , 1 001  "  v 2c ^ n t ) 


pqst, 1010 


] 


(1.1-19) 


and 


% 

L. 


in  which  c -nt  is  the  internal  heat  capacity  per  molecule.  The 
spherical  components  of  A  and  n  are  then  given  by 


The  important  point  of  Hunter's  development  of  the  transport 

coefficients  is  that  the  tensor  analysis  precedes  any  truncation 

of  the  basis.  Solutions  to  (1.1-19)  and  (1.1-20)  require  the 

truncation  of  the  basis  and  evaluation  of  the  scalars, 

o  S i ^ i ) ^  •  The  remainder  of  Part  1  is  concerned  with  the 
p  q  s  t  j 

evaluation  of  these  scalars.  With  the  ultimate  goal  of  obtaining 
calculated  values  of  these  scalars,  trey  are  expressed  in  terms  of 
relative  momentum  collision  integrals  in  section  1.2,  which  in 
turn  are  written  in  terms  of  the  reduced  scattering  matrix  of 
Curtiss  and  co-workers^  in  sections  1.3  and  1.5  and  simplified 


by  the  introduction  of  various  operators  in  sections  1.4  and  1.6. 
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Finally,  some  interesting  relations  among  the  scalars,  suggested 

4 

by  similar  relations  found  by  Chen,  Moraal ,  and  Snider  for  the 
collision  integrals  in  another  representation ,  are  developed  in 


section  1.7. 
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1.2  Expansion  of  the  0 |^p •  q ■  5 •  t * j 
Momentum  Collision  Integrals 
(a) 


(a) 


in  Terms  of  the  Relative 


The  o 


p  q  s  t 
p'q's't1 


provide  a  description  of  the  collision 

process  in  terms  of  initial  and  final  momenta  (psp's1)  of  a  single 

molecule.  In  order  to  express  these  quantities  in  terms  of 

S(j‘j‘A|S  S.  )  ,  the  reduced  scattering  matrix,  it  is  first  necessary 
a  d  a  d 

to  make  the  transition  to  a  description  of  the  collision  process  in 

terms  of  initial  and  final  relative  momenta  (xn^-'n1)  of  the 

collision  pair.  This  is  accomplished  through  the  introduction  of 
4  ( k  1 

the  quantity  I  ,  ,  ,  ,  ,  defined  in  CMS-B16  and  resulting 

*'1**'  n  ) p s p  s 

from  the  transformation  to  center  of  mass  and  relative  coordinates. 


This  section  is  concerned  with  obtaining  expressions  for  the 

-IS-SiM**'  <"  te«  of  th*  ’'(c'qvf')k  and  0”('v!!'r)k 

relative  momentum  collision  integrals. 

The  scalars  a j ^ ^ ^ ^ ^  defined  in  H-40  and  (1.1-15)  are 
related1  to  the  scalars  r(pqsi  (p'q's1 j' ) ^  defined  in  H-24  and 

a  a 

(1.1-13): 


ofP,  =  l  «“1(qq*  )(p-  ([.I](q)'-:9[J](q))^ 

lP  q  1  J  j  j’  JaJa  ‘  ' 

3  d 


<([-!](q,)^,[J](q,))’/(  R[q>(>-i  )R[?,)(f-i.)  o(pqsj  a|p'q's*  j*  ){a) 

Ja  Ja  Ja  a 

(1.2-1) 


i n  wh i ch  a  ( k ;  <  :.  *  *  ’ )  ~  *'(2  ■  j4 1 )  (2C2+1 )  •  •  •  , 

C:  =  (E.  /kpT)  =  (fi2j  (j  +1  ))/(2I  k„T)  is  the  dimensionless 
j  j,  o  a  a  a  0 

a  d 


rotational  energy  of  molecule  a  ,  and 

(p,  i.)I/2  =  (P.  )  /?(P;.)^  =  ^ ( J, j '  )Q_1  exp(-e  -  -cj.)  ,  Q  being 
JaJa  Ja  Ja  Ja  Ja 

the  molecular  partition  function.  The  quantity 

([J](q)^]{qV/2  is  introduced  in  Ref.  4. 

Ja  la) 

Using  H-27  and  HS-19,  the  scalars  ofpqsj^  Jp'q's' j^)  '  can 

be  expressed  in  terms  of  the  scalars  o(pqs j  |p‘ a’ s '  )(<  ,  which 

in  turn  can  be  expressed  in  terms  of  the  cross  sections, 

a,(«njajbq|e,n,j;^q‘)k  and  a"(  ^jaJbq  |  t'n*  j^q1  )R  . 

Substitution  of  these  expressions  in  (1.2-1)  yields 


a  P  q  s  1 
p '  q '  s  '  t  * 


1(a) 


(-l)q+d+P'  l  (-1  )kct2(k)n(kq*q)  y2^kp.p)  ^2r,(kpp* )'  ^{2-  p‘  a} 


”  L  4  "Wr* 

2  >  n  1  a  D 

x,  n  i  >  ,  > 

Vb 


(CJ](q)^CJ](q))!/2  ([•I3(q,)d1'[J](q,))'/2 


X  R[q)(e,  )  r{?’  ){ci. )  a’^qq1) 
Ja  Ja 


{o,(tnjajbq|«.'n'j;j'q,)k  +  (-1)£  a"  Unja  Jbq  ( v.'  n 1  j  ’  j^q* 


(1.2-2) 


The  quantities  n(K.18.2?'3)  are  defined  in  CMS-A22  and 
{p1  p’  a}  is  a  6~j  symbo1-8 


H-63  defines  a  new  scalar  which  can  be  written 


X 


a^Mkp'pr'M 

£'n' 

l  ( P i  i  i,i,)V?([J](q)^[J](q))/2  ([J](q,)€P,[JJ(q,))^ 
j  j.  JaJbJaJb  Ja  Ja 

3  D 


X  R [q)(Cj  )  r[?'  )(cj. )  cT^qq1) 

3  3 

x  {a- («lnjajbqK‘n*  )R  +  {-])V 

This  can  be  rewritten  as 


c’^tnj^qlrn'j’j-q')^ 

(1.2-3) 


p  q  s  t 

p '  q ’ s  ‘  t* 


v(k)«kp'pr‘^  £j  «u<>«  i^.„.;psp.s, 

t’n’ 


*  «cj 


P  q  n  t 
rq'n't'J 


+  (_l)e'  a"|£  d  n  t  | 

'  '  It'q'n't 


}  (..2- 


4) 


if  two  new  collision  integrals  are  defined: 
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.  U  q  n  t 


*‘q‘ 


n’t*  j. 


•  “  ■  lPj  j  i * i ’ ^ 

a  JaJbJaJb 

a  D 

jajb 


ly2 


(CJ](qM[j])(q))!^  (CJ](q')^‘[3](q,))|/? 


x  cf^qq'jR^Uj  )r£?  ^(ej.)  a’  ( ZnJa  jfaq  1 2,'n 1  jjj^q'  )k  (1.2-5) 

3  3 


and 


4*  q n  *  1  =  <•  (D  ) 
i,qVt,Jk  ’  XWt 

a  b 
jajb 


V, 


j  a  ~  -  J  a 

a’WjR^U.j  jR^’^ej.)  a" ( Znja jbq | Z* n *  )k.  (1.2-6) 


These  are  the  desired  relative  momentum  collision  integrals.  The 


>  q  s  t 
[p'q's't* 


pqst 

p'q's’t1 


(a) 


(a) 


scalars  can  now  be  written 


l  n(kq'q)lA-  r.(kp‘p) 1/2  (-l)k+q'a+P' 


q  q’  k]  jp  q  S  t 
P'  P  (p’q's't'J 


=  (-l)q+a+P‘  y  (-1  )k  ^(k)  r.(kq-q)1^  |q,  q'  k 


ipsp's' 


*  |n  S*k«')'4  tini-n-,, 

l'n' 

f  ,  ft  q  n  t  )  ,  /  ,  ,,'i  q  n  t  1  \  ,  T  7S 

“  i’  h'.'t'li,  ‘-1’  °  U'q'-'fjJ  ■(L‘-7) 


£  q  n  t  1 

After  lengthy  substitutions,  the  o'  ^ • q • n • ^ ^  ar>d 

relative  momentum  collision  integrals  can  be 
expressed  in  terms  of  the  reduced  scattering  matrix,  S( lsaSb^ ’ 
and  consequently,  in  terms  of  the  generalized  phase  shift, 

H(  j  1  A I  S  S.  )  ,  of  Curtiss  and  co-workers.^  Considerable  algebraic 

d  D  3D 

manipulation  and  the  introduction  of  various  operators  allows  the 

expressions  for  °  I  t'q*  n*  t*  Jk  and  0(rq,n,t‘lk  t0  be  greatly 
simplified.  These  substitutions  and  simplifications  are  carried 
out  in  the  next  four  sections. 


.,  [£  q  n  t 
t'q'n't 


1 

l 
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1.3  [j'q'n1  t‘)k  'n  7erms  of  the  Reducet*  Scattering  Matrix 


Equation  (1.2-5)  is  an  expression  for  a' 


q  n  t  ] 
'q’n’t'Jl 


in 


terms  of  the  kinetic  theory  cross  section,  a'  ( 2.njaJbq  1 1‘  n'  j^j^q*  )k 
for  nonvibrating  diatomic  molecules.  HS-43  gives  an  expression 
for  a,(enjajbqil'n'j^j^q,)k  in  terms  of  the  relative  velocity 
cross  sections,  oi(£n| j a j ^ ; kO j ?'n' )k  and 

^(tnjAlWb\;LaNaLbMb'1'n'jajb)k  *  s“bstitution  °f  this 
expression  in  (1.2-5)  yields 


iqnt] 


=  f 


:(kqq')~1/2  iq+q‘ 


[  (  Dk  o‘1(k)a(j-)(p.1.1) 

■jajb  JaJb 


«. 1  q '  n '  t  ‘  J 

-  ([s?](q)^[J](q)).y?  ([J)(q,)©q’[J](q,))^  R[q)(e  ,)R[q,)(c  ,) 

J  a  J  a  L  J  a  L  J  a 


3  3  i[o  j  ( P.n  |  j  ’  j  ‘  ;k0 1 2.*  n* ).  +  (-l)q+q  o(  In  |  j '  j  ’  ;k0(  4*  n’  )*  ] 

q'  j!  k  a  0  k  a  b 

a 


l  *}(ihU2(yMiLHp<  , 


V, 


■Vb 


K  '  '  '  i  i  l*!1 

°  JaJbJaJb 


([J](q)©q[J](q))^  ([.]](q,)©q,[J](q,))^  R[q){t:.  )r[?')(c.,) 

A  J  a  ^  ^  J  a 


LaLaLb 

M  M’M. 
a  a  b 


f<J*  q  k  } 


^a  ^a  La 
^a  ^a  La^ 


a(L  i:)(-l)La+Mb 

a  a 


La  l; 


a  Ja 

k 


M,  -M'  M*  -ft 
a  a  a  a 


"J(«njaj1)|LaMaLbH;L;HiLbMb|i'n'Jiji)k 


U. 3-1) 


q’  q  k 

in  which  ij‘  T  l  •  is  a  9-j  symbol.1 

add 

i'  i  L' 
xJa  Ja  aJ 


The  relative  velocity  cross  sections  are  expressed  in  terms 
of  quantities  F, (j' j* ;L  H  L. M, ;j  j.  )  with  the  aid  of  HS-41  and 

A  d  D  a  a  D  0  a  D 

HS-42: 


oi(tn|j^j^;kO|e'n‘ )k  “  \  j8J^T 
x  fk  f  *')  f  hv  yl .  .t 


n(k£f  )_1/2  (-i)S,+i'  (-Dk  a(U') 


o  o  o  j  J  dY  l  WA>k00°->W 

(1.3-2) 


^^aJbiLaMaLbHbsWbHbl1,n,JiJU  =  IT  )_ 

x  [  dyy2  [  dg‘ (g’ )2  ;?  (y)r  (y1  ) 

M  -M‘  0  M'-M  j  99  nl  "  1 

a  a  a  a-' 

fx  X'  £' )  f  X  X'  £'  ' 

x  I  ct2(XX’)!  j 

XX1  10  0  C  M  +M.  -M'-M.  M'-M 

v  v  a  b  a  b  a  aJ 


* ,Wb;WbVVb>  ■  (I-3-3> 


In  the  above,  g'  and  g  are  the  relative  velocities  of  the 
centers  of  mass  of  the  colliding  molecules  before  and  after 
collision,  while  y'  and  y  are  the  dimensionless  relative 
velocities  (momenta)  before  and  after  collision.  The  R^y)  are 
normalized  three-dimensional  harmonic  oscillator  wavef unctions 


introduced  in  CMS- 74. 


The  Fj  quantitif 
sections  are  normalize^ 
potential ,  V  =  , 


; pearing  in  the  relative  velocity  cross 
ich  that  for  a  spherical  interaction 


^A^a^b’^a^a^b^b  ’^a^b^spheri  cal 


<  [1  -  expUir,,)]  ,  (1.3-4) 


nx  being  the  spherical  phase  shift.  Thus,  the  are  a  type 
of  reduced  transition  matrix.  In  terms  of  S(j' ji A  IS  S.)  ,  the 

3  D  a  D 
5 

reduced  scattering  matrix. 


aJb; 


“a  b1' 


x  (8tt2)~2  l  (-1 

jV+p 

pa  ^  Lal 

4  -^b  Lb 

vp 

"P^ 

1 

P 

o 

0  u  -U, 

JJ 


dSa  dSb  [,  -  S(JXMSaSb)]D^  (Sa)  oh*  (S„)  .  (1.3-5) 


where 


1  L  1 
a  a 

0  v  -v 


i  i 1  L 
Ja  Ja  a 


is  a  3-j  symbol  and  the  (S)  are  elements 

vM,  3 

8  a 


of  the  rotation  matrix.' 

Insertion  of  equation  (1.3-5)  in  equations  (1.3-2)  and  (1.3-3) 


leads  to  expressions  for  the  relative  velocity  cross  sections  in 
terms  of  the  reduced  scattering  matrix: 
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O,(tn|j;^;k0|£’n')k  =  fvL(s^r)/;  I“2(X) 


[fXO££)"1/2  a (£)  <S(£k|£'0) 


-  [n(m*)'1/2  (-i)£+r  {-l)k+ja  a(j;k)a(££')[o  J  J‘ 

dVwRra(yK'V(r)  l  {~i)V  (o3  va  -v)  (8jr2)'2 


dsa  dsb  S(j;j'x|sasb)  Dk0(sa)] 


(1.3-6) 


and 


fi2ir 

ITU 

U 

[8kBTj 

0  2  ( xn  J  a  J  b  I  LaMa1'  bMb  ’  LaMa^bMb  ^ £  * n  * J  a  J  b  ^  k  = 

x  |  n(kZi)'^  (-1)Z  a2(£)[k  J  *] 

-  6(£i£,)6(jajb|j;j')6(LaL;Lb|000)5(MaM;Mb|000) 
x  I  a2(X)  !  dY  ^Rn£(y)  Rn-£(Y)] 

-  jr.(k££‘ (-i)£+r  (-l)ja  a(U’)o(j;)  6(  JaJb  |  1 00 ) 


|dY  ug  Rn£(v)  Rn’£,(Y)  a^AA^l0  0  0  ^ 

A  A  V 


(8tt2) 


-2 


■  r 


dS,  dS. 
a  b 


"a  a 


VU  c(Lj 

X  a' 

£*' 

k  £  £.' ' 

(j1  y  l  ) 

Ja  Ja  a 

a 

M  0 
a 

-M 

aJ 

M  0  -M, 

'■a  aJ 

0  v  -v, 

s  IWW  o^J  (s4)  )]*Mi/a|oo)H)'VV;> 


[X  X'  £ ') 

M'  0  -M' 
a  a 


ik  £  rU y  y  lm  * 

add 


M*  0  -M' 
a  a 


!  0 


_>  ' 


!  s'iWis.VDj',s.:‘ 

d 


(Eauation  continued  on  following  page.) 
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jaOarrVi)i+r  a  (IV)  |  d>y2  f  dg‘(g‘)2  R^y)^,  v  h' 


X  V  l '] 

(  X  X1  V 

k  i  V 

0  0  0  J 

) 

[m  +M.  -f1'-M.  M'-Ma 
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M  -M*  0  M'-M 

a  a  a  a-* 

*  l  a2(XX' ) 

XX' 
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dSadSbdS;dSb  S<J^bXlSaSb>  5*(Wls;s^) 

(1.3-7) 


-  <v  °jV<sb»  D'Vs;>*  °uV(sb>* 

a  b  a  b 


At  tbs  stage  it  is  useful  to  consider  the  terms  of 

iLV,^,  involving  relative  velocity  cross  sections  as 
[t  q  n  t  )  ^ 

they  appear  in  equation  (1.3-1): 


fUkqq1)'^  (i)q+q‘  l  (-Dk  a'1(k)c(j;)(p.1,1)([.j](q)9q[j](q))/; 
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(87i2)'2|| 


}■ 


since 


8,9 


dSadSb 

(1.3-8) 


q?  j.  0a|  =  (-Dq  a_1(j^q)  6(q|q')  .  The  and 

4 


summations  appearing  in  the  first  term  above  can  be  carried  out: 


l  (Pi.i,)([J](q)^1[J](q))i,  R|q)(ei.)R[q)(c.,)  =  u2(q)  6..,  , 

;i  ii  “  Ja  J-  L  J, 


j']'  JaJb 
JaJb 


(1.3-9) 


allowing  the  right  side  oi  equation  (1.3-8)  to  be  written  as 


=  ^(g^T]  ■  I  u2  ( A )  |  [ n( Oqq )  ‘  !/;’  fi(  0£ £ ) " 1/2 
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*  (P  4  .  t  .  ) 
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V? 


o  0  0  j|q'  j*  k  j 
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(1.3-10) 


The  above  expression  involves  a  constant  term  and  a  term  linear 
in  the  reduced  scattering  matrix,  a  result  of  the  fact  that  the 

5 

cross  section  is  linear  in  the  transtion  matrix. 

The  o2  cross  section,  on  the  other  hand,  is  bilinear  in  the 

5 

transition  matrix,  and  consequently,  this  term  is  somewhat  more 
difficult  to  evaluate: 


fi(kqq'f  V2(i)q+q,H)q  [  a'1(jb)a2(j')a(j')(pi  .  . ,  . ,  ) 
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(-i) 


(8,2) 
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I 
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> 

1 

O 

[o 

u  -Uj 

(0  v' 
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[dg'(g‘)2S]-^,(Y)R„.,.(Y,)(8^)“+  f[jf  ds,dskds:ds^ 
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L. 


a  b  a  b 
L 1  L 


«  s(j;j;Ms,sb>s  u^is^io.;  <s,>ox<volS.<s;>oA<s> 

a  d  a  b 


(1.3-11) 


In  the  above  expression. 


k  i  l 

0  0  0 


q 


o  0  kl 

0  0  0 


(-1  )z  Ct’1  (?.q)cx'2(j'  )  5 ( q k  ; q ■  0 )  , 

a 

(1.3-12) 


allowing  the  ja>  jb>  j^,  and  sums  to  be  carried  out  in  the 

first  term  in  brackets.  The  A'  sum  can  be  evaluated  in  the  two 

8  9 

terms  in  the  second  set  of  brackets,  ’  yielding  o(Maj0)  and 

o(rr  [o)  . 
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fq‘ 

i 

q  k 

(L  Ok] 

j  j  ■ 

j 1  l 

a 
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M  0  -M 
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t  a  a-1 

(-1) 
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d 


y'\ 

Ja  j 


q1  ji  k 


(I. 3-1 3a) 


and 


q  q 

,ja  ja 
ja  ja 


0 


k 
0 

LJl° 


La  k 
‘Ma 


(-l)q,'Ma  a'2(L;)r'(j-) 


fj]  q  j 

6(LMk)  a 

W  Ja  k  i 


(I. 3-1 3b) 


allowing  the  evaluation  of  the  j,,  j.  ,  L  .  M  ,  L’,  and  M'  sums 

d  D  d  d  d  d 

in  the  terms  in  the  second  set  of  brackets.  The  right  side  of 
equation  (1.3-11)  then  becomes 
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(1.3-14) 
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Considering  these  results. 
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★ 


(1.3-15) 


Before  proceeding  to  the  reduction  of  this  collision 


integral,  it  is  useful  to  integrate  over  the  energy  delta  function 
appearing  in  the  term  quadratic  in  i  .  When  this  is  done. 
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(1.3-16) 


Equation  (1.3-15)  can  now  be  written 
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(1.3-17) 


an  expression  that  provides  a  convenient  starting  point  for  the 
general  reduction  to  be  accomplished  in  the  next  section. 


r*<rvf  ~r.'  fWh-jfcra  n’  *7 
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1.4  Reduction  or  the  o'  - , ^ ^ ^ ^  Collision  Integrals 


In  this  section,  equation  (1.3-17)  is  reduced  to  a  summation 
over  eight  indices  and  an  integration  over  six  angles,  in  contrast 
to  its  present  form  requiring  a  sixteen-fold  summation  and 
twelve-fold  angle  integration.  The  reduction  is  made  possible 
by  the  introduction  of  certain  operators,  as  well  as  the  explicit 
summation  over  several  indices. 

In  Reference  10,  it  is  shown  that 
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and  as  a  consequence, 
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,(q) 


Remembering  that  R|h'(‘  .  )  is  a  Hang  Chang-Uhlenbeck  polynomial. 


it  is  possible  to  write 
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where  the  operator  is  defined  as 
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0  V 


'a  “a  'a  II  Jb  Jb  “b  J~a  , ,  ,  J  b 


-vj  io 


1  "vK  >V  lV 

.J  -uj  a  h 


I  (-1) 

vu 


v+u 


4  4  La 

0  o  -v 


4  4  4 

0  u  -vjJ 


J  dy  y‘ 


*’  e-'-  Ls:*l/j(v2) 


dSadSb  !**»'«•*>  S<j;ii*|S,Sb))  oJJ  (Sa)  D  d  (Sb)  . 


(1.4-3) 


This  allows  equation  (1.3-17)  to  be  written 


,[t  q  n  t  |  frjn 
'  ( k”  q 1  n  *  t  ’  j  ^ 


4>jkJ-  {[[  «?(A)  J2(0qq V'h  k.(0..)'V-’ 


c(vq) 


v  n( 1 qtk  J  Z *  q *  t 1 0 )  f  d> ,  »„,(>)«„. f (y)J 

2i;(kqq,)'l/-  !;(kLf  )"'/;'  ( - 1  )q‘q'  4  e+£'  a[U' )  l  £  J 

44  LaLaLb  AA 
44  MaMA 

<  .  >Ma‘Ma  ,  .  iv+u-v'-u'  ,  ,>l-a+Mb  3,  .  .  ,  .  .  , 

(-1  '  (*’ )  (-1 )  (i J(Ja)  a-  (jb)  a(j^) 


(Equation  continued  on  following  page.) 


-  {pi,i.)([J](q)^[J](q))iA  ([j](q,)^‘[j](q,))/’  r(L  L'L.  ) 
-Vh  ~  '  -  J,  aab 


fA  A'  V 


0  0  0  JIM+M.  -M'-M.  M'-M  I  M  -M*  0  M'-M 

'  ^  3  D  3D  3  3'  ^  3  d  3 


k  i 


[q 1  q  k  1 


1  v  i  |  1 

tJa  Ja  La| 


4  iJi  b  L;i 


La  La  k  ]  ja  ja  LaMja  ja  La  jb  Jb  Lbj  [Jb  Jb  Lb 
[\  -M'  MM>lj[o  v  -«,j[o  V  -o'  0  i  j  |o  p'  -u' 


:(nH  )  P  ( n 1  + 1 ) 

F(n+?  +  ^2Tr  ( n '  + 11. 1  +  %  ] 


'"T4  f  dw‘  e'r  L^.  f‘/?(r  )R^q  Hey) 


dSadSbd5;d^  £Cab('qnt>  ^bAlSaSb>l 


1-3  K  1  K  \  nLa  KM*  nkb  KM*!  . 


xS  is;sb>  °vSa(sa)  °;.tVsb)  °vV,(s;,)  D\^j] 


(1.4-4) 


Now,  the  sum  can  be  evaluated/ 


l 


|Jb  jb  Lb] !Jb  Jb  Lb  1 


[0  u  - 1 1  j  i.  0 


=  6  ,  , 
uu 


(1.4-5) 


and  the  9-j  symbol  expanded. 


M  n  t  I  a  JILL-Jf  ^(x)-2(0qq)"  ':'S2(0^.)  «»  Aq) 

rq,n’t,]k  4uk0T  Vv 

i. 

X  6(£qtkl a* q’ t’O)  [  &,  Y^n?(T)  Rn ■ e 1  > 

_  1/  _  1/  q-q 1  t  i-1.’  ‘  , ,  , ,  v  j  ,  'i  *  1  J_  1  (jilt)  J  I 

-  zn(kqq‘)  ^(kU*  )  4  ‘)i  Ylnn,^.,  +  7)}i 


(3tt'T  >  l  (-») 

i  I'i!  L.M  M  H'  -  1 1  j 

JaJaJb  b  d  a  a  .  {;_•  , 

•  vyv' 

U\  1^2 


V'V'1 


<+3a^a+“i+^r+Mi"+Maf% 


d  a  J-Ju 


J  A 


A  A1  £'U  A 


k  i 


o.UvMb  -M;-Mb  M;-Mai'^d'M;  0  M*-Ma 


■fja  K 

Ja" 

,  lo  1 

q'  Ja  jaMja  q  ja  \  \  K  ja  ja  3  q  | 

Ui  v’  V2.  (0  U2*  -M;' 1  .  .-Ill'  v  >>J  2  ,  1  'iJ2  ^a  ^ai 


’  yfl  e-f2  l^,^r’)R<?'>(Ej.)JJ||dSadS|)dS'<i5' 
[C*b(!lqnt)  su;j’»|sasb)]  S*( j;jjv '  1  s;s'  )D^  ( y  °^<y* 


(Equation  continued  on  following  page.) 


a 


L ' 


U; 


~V2  ' 

-j 

l(Ma  u 

r" 

(q1 

)| 

K 

q  l; 

i-M 

-  VI  - 

.  -v’ 

-Ma 

i  D  (S  ) 
h  j  oM  a 

a 


LI  * 

D  a  /S') 
*VM'  Ua; 


(1.4-7) 


Using  Ref.  8,  the  sums  over  L  and  L'  can  be  carried  out, 

a  d 


’  K 

q 

L 

a 

La  q 

M\  * 

~iJ  2  * 

-V 

Ma  *Ma"'1?"- 

L. 


D  “  (S  ) 
vM  '  a' 

a 


v+M, 

(-1)  3 


k  q 

‘'v.m'-ii*'  Duj',Ma+u?,,(Sa)  ,-u2"(Sa)  ’ 

(I.4-8a) 


and 


L1 


L '  11 

l  (-1)  fl 


a 


K  L 1 

'  q’ 

L1  r  ) 

a 

a 

-Hn  -v*. 

.-m" 

-M’  M  +uo" 
a  a  - 

dv'm'<s;»’ 


l’+  r 


(-1  )q  '  1  .  6.,,  M  „  Dq  AS')  D  M  n(S 

v  Mfl  » Ma  ~  M  ]  +  M?  a  *‘2  • 


(I.4-8b) 


followed  by  the  sums  over  v,  v’,  and  Mr  .  After  relabelling 

d 


indices , 
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,  fi  q  n  t 
Z'q'n' t*  J 


"h2:r 


A 


/, 


4Tk-  T'  •  [/  a2(A)r.(0qq)'  2;'.((U;t)’ '  i(J.q) 


B  i  A 

l 


x  5  ( Hq  t  k  |  S.  *  q '  t'O)  dy  y  Rn£(y)  ^(y)] 


-  2P.(kqqT^  tfkUT*  »(n-)i^ii;tJian&H 


i 


A 


i‘\n+i+->5)r(n,+s.'  *}{Y\ 


(8tt2)'4  V  !  AO''1  )a2(L--) A(  j  )nt(j')(p 

V^b  ■-  *  ’  ,.jJb 


'£ 


(cj](q)^[j](q)i.'  ([jj(q,)^,cj](q,))i; 

4  U  j 


VJ 

I 

a 


l 


(.i)Mi-ur-u2+P2,+b3,,-M2"  (.1^a+ja^2+yi,+M1‘ 


Ul,U2,Ml"U2,,U3" 


X 

X' 

V ' 

f  X  X' 

V  ' 

'i  k 

v  1 

0 

V. 

0 

0  _ 

1 

1  .  ,  .  l  II 

rp  u+;.2  -  u  3 

3  "■J2 

0 

u3"-P:" 

'ja 

q 

|  *  j; 

4' 

■Ja  < 

j;  i 

j 

o' 

j; 

Ja  1 

o 

U2 

-U2, 

I'bj  Uj-u 2 

U  ;  _ 

1°  U2* 

•u,1  j 

i  ’ 

.  i  ,  i 

Ul  “42 

u? '  j 

k  q  q' 
.U3"-U2"  U2"  ~M  3 " 


I  dvrUI  e'y; 


*  Jj  (S,)0J  „  .(S,) 


(Equation  continued  on  following  page.) 
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X  Drb  *  "  ,  «(SJ  i  I 

S.-p+U]  ~  VJ;»  D  jj 


m  s  'W'lW 


^  Dqt'  ,  ..(S' )  oLb  ..  ,.(s;)’ 

-  U  i  -Mj  a  -(j?  D 


(1.4-9) 


At  this  point  it  is  convenient  to  introduce  a  compact 
notation  for  the  i  summation; 

a 


Kq^jKjiyiw?')  =  l  (-DJd  <v’(ja)([J](q)^[j]{q))./? 

ja 


(1.4-1 Oa ) 


Ja  q 


(0  Vi  2  -U2j 


K  J. 


J  a  I  1  K 


a  “a | |  a  | 

-Vil  i1 1 ~ l-1 2  i|;>j(p2'  “I1?'  0  J 


as  well  as  two  operators  defined  on  the  S  S.  and  S'  angle  spaces, 

(3D  Q 

respectively,  by  the  following  eigenvalue  relations: 


(Sa)D“  (Sb)  .  (-1) 

d  D 


-FT ’.Mi 


!  j  D^.  (S  )  D^.  (S.  ) 

|fT%  -f,M.M3+P-Ma  Mfl-r!  ’’'a  a  'iMb  b 


(-D 


M,-M.  M.+f.  M  -f: 

3  b  u  d 


“vM  Isa)  DX(S»> 

3  b 


( 1 . 4-1  Ob) 


and 
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La 


3 


H) 


(a'} 


h~  1  L 


a 

(a1) 


"i‘t 


a 

-h"  1  L 


i  /  (a * ) 


DLa  (S') 
v'M’ v  V 
a 


(-D 


j.>’fq'  K  J,!  L 


U;'  v'-Ui’ 


n  (s’  ) 
U/M'1  a'  ’ 
a 


(1.4-1 Gc ) 


(a 1  ) 

where  M3  and  L\  ’  are  both  defined  in  Ref.  10.  Now 


Lb 


E  b(^,i,;u3"-u.»")[D  „  „(s  )Dm  „(S  )D  u  „  „(S  )] 

ao  m  1 —  M 1  a  'II2U'!  a  ii,-p+ui  -u;>  b 


=  (-D1 


X  V  V 

-M  U^2h-M3"  U3"“U2 ' 


[D  *  „(S  )Dq  ,i  ( 5  ) 

u  1  -  u  1  a  -  j  1  u  L  a 


Lb 


X  DB,-u+ui"-ii;'"(Sb^  ■ 


(1.4-11) 


If  the  3-j  symbols  on  the  RHS  of  equations  ( 1 . 4- 1  Ob )  and  (1.4-1 0c ) 
8  9 

are  expanded,  ’  they,  along  with  the  phases,  become  polynomials 

in  (Ma+Mtj)  and  v1  ,  respectively.  Similarly,  the  operators 

(a 1  ) 

and  E^,  can  be  expanded  as  polynomials  in  ,M3  and  .  It 

should  be  noted  that  K  and  ,M3  commute,  while  K  and  ' 
do  not.  Thus  K  commutes  with  Eflb( XX1 £’ ;p3"-p?" )  but  not  with 
Ea.(q'j’;pr)  . 
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Once  the  value  of  q  has  been  specified,  the  I (qp2 ) 

a 

sums  can  be  evaluated  using  recursion  relations  among  the  3-j 
coefficients.  Values  for  I (qu? j Ip iW? ' )  are  given  in 

a 

Appendix  l.A  for  q=  0,1  ,  and  2.  These  are  sufficient  to  evaluate 

transport  properties  in  an  applied  field. ^ 

Introducing  the  above  and  making  use  of  the  fact  that  M3 
(a 1  1 

and  L'  '  ,  and  consequently  Eajj(AA' 2.' ;p3"-p2")  and 
£  ,  (q1  j  ‘  ;u  i  * )  »  are  hermitian,^  it  is  possible  to  write 

d  d 


M  n  t  1 

fi2^ 

2,q'n't'J 

k  4ukBT 

U  ct2(A)fi(0qq) 

A 


2ft(022)  '2  a(2q) 


V2 


x  6(2qtk|2'q't‘0) 

-  2$i(kqq’  ) 


dY  Y  Rnt(Y)  Rn.t(Y)] 


-K2(  1.)q-q,+2+2'  h(££() 


(8~‘)‘4  y  [^(AA'j^d hK)«(v)(P i<ii)([j](q,)^,[j](q,))f? 

j; Jb  b  "  JaJb  ‘  "  Ja 


XV 

V 


l  (_i  )bl_Pl  '  -U2+Ul'  1  +M3"-M2"(_1  ,+Ml  ‘  +M  1  ‘ 

PpUll'C 


U)  ’p?'p]"p?"p3" 


X  >'  2'1 

'  •  k  2  ’ 

{  k  q  q' 

0  0  0  , 

(0  pV'-ub"  p3"-p.". 

,  II  ,  II  .  .  M  ,ii 

(U3  -V?  U2  -p3 

(Equation  continued  on  following  page.) 
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I(qM2;KjaUn,2,)R^  j  dy  v’+1  e’Y<>  L^/  (y2) 

{jf  «,«bte*(>v*-!„,--M.*)4(w)stj;^i5,sb)]iiJ|.  .<s.) 


Dq  „(SJ  0gb  (5.  )i 
-M2V1?  a  3,-;.  b'J 


*  “u/u,-'5^  °tu<sb>*[ 


M.4-IZ) 


In  order  to  proceed  beyond  the  above  expression,  it  is  useful 
to  introduce  another  operator,  I  ,(qp2)  ,  defined  by  the  relation 

a 


i>m2>  f(s;)  =  (8^)-'  l  0?  ,.(s;>*  0)'“^' 

<li  j 

Pi'V  { 1 .4-1 3a) 


x  I(qP-;Kj ■ui^2' )  f  dS"  D<  ,.(S")  f(S')  , 
a  j  a  a?  ji  a  a 


where  f(S')  is  an  arbitrary  function  in  the  3  anyle  space. 

a  3 

As  a  consequence  of  (1.4-1 3a), 


l  _  (i)Ul+u"'  j  ds; 

V2 

-  .  dS;  f<s;)  *  f  «5;  f(Sa)la,(qll;)D.;iiii..(S;). 


( 1 . 4  - 1 3b ) 
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•j-  * 

taking  the  adjoint  of  I  , (qp2 )  •  The  above  expression  can  be 

u 

used  to  determine  I  ,(qp2)  Tor  a  given  I(qp2;K j'P]U2)  .  The 

a  a 

Za,(qp2)  corresponding  to  the  q  =  0,1,  and  2  values  of 
I(qp2-,K  J  'piW-3  1 )  ace  also  given  in  Appendix  I. A. 

cl 

The  7  . ( qu2 )  operators  isolate  all  the  explicit  dependence 
on  the  indices  and  p3  appearing  in  the  I ( qv^ 2 J 1 M 2 ‘ ) 

K 

in  the  rotation  matrix  elements,  D  „(S‘)  ,  instead. 

-PiPl  a 

Equation  (1.4-12)  can  now  be  written  as 


17  q  n  t 


Q  a2(X)Q(0qq)"1/^f40£J:)'1/2a(i?.q) 
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K,q,n,t'Jk  "  4ukgT 
X  6 (> qtk | £' q’ t'O)  J  dy  y  ^(y)  Rn«  £(y)3 

-  2S2(kqq,)'^r.(kU')'^(i)q"q'+f+r,'t2(£t') 


.  T(n+1; 

r(n'+l) 

rrii+I+H 

r(n'+J7,  +  ^) 
— - 

(8:-2) 


-4 


Jajb 

A  A  * 


1  \  H 


o?(>.A,)x(j;)(p.,,,)([j](q,)^,[j]{q,)); 


a# 


l  (_i  )'-‘3  ,+V‘2f‘ 


Up  II] 

u"  1." 


X  1 

£' ' 

£  k  V 

0  0 

0  . 

0  Mr."-P,"  P3"-P2"; 

II  II 


1 1' 3  _P?  U2  _P3 


q  q'  I  \ 

’(c,.)  I 

"  1  J 


a  X.+1  __Y2  ,  +  1/, ,  7  < 

dy  y  e  L  ,  ‘(y2) 


(Equation  continued  on  following  page.) 


(Equation  continued  on  following  page.) 
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Further  simplification  of  the  above  expression  is  difficult 
without  restriction  of  one  or  more  of  the  indices. 
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1.5  o" 


IWIU  in  Terms  of  the  Reduced  Scattering  Matrix 
a.  q  n  t  i  K 


Equation  (1.2-6)  is  an  expression  for  0 "  ]  v  ■  i j 1  ij !  t ’  i k 


'  n 


terms  of  the  kinetic  theory  cross  section,  u"(  -  nj  i |  1  n‘  j^q ‘  )k 
for  nonvibrating  diatomic  molecules.  HS-44  is  expression  for 
o"(  j  i‘  n‘  j.’j  ( g '  in  terms  of  the  relaLive  velocity  cross 

sections,  o  x  ( j  jjs)b;qq  '  ’  •  ‘  n*  )k  and 

!  .(tnj  j.  |L  M  L. :<t.  ;L  'M’LuM'  j  «  ' n'  : '  j  '  ?  ,  which  ;m  Oh  n  are 

■v  JaJb'  a  a  b  b’  a  a  b  b1  -b  a  r. 

expressed  in  terms  of  quantities,  E.  ,  using  HS-'ii  and  HS-4?: 


di(«nl  jajb;qq'  I  ?-'n'  )k  =  t(-o'k) 


k  1  i  • ) 

0  0  0  J 


;  dv  ^  R  p ( v )  (y) 

)  pg  nx  fi  <. 


q'  k) 


AM 


l  u:'(A)  |  F  (j  :  ;q  M  q'  -M  ;.j  t  .  ,  (1.5-1) 

,,  r, ;  \  a  n  a  a  a.  .. 


If1  -M  Oj 
a  a  • 


and 


o,(lnj  j.  |L  M  L.M.  a’M'L'M;  I  •  '  n '  ij  j'V 
1 v  Ja  b 1  a  a  b  b  a  a  b  b 1  b  a  - k 


fr’  rr  f  - 
p  j  Ht  .  i 


-■(fc  ■  -V  ) 


-  V 


<  (-i)5><'  X  ( ‘  )  |  d  y  y  •  j  dq  *  k  g  ‘  i  Rn.(,  •  Rn,  .  .  (l‘  ) 

f>  a'  i'H  '  A* 

X  l  v’(AX,)l  !■  j 

AA 1  in  0  0  I  |M  ifl.  -M'-M’  M'+M'-M  -  M.  I 

'  '  a  b  a  b  a  b  a  b' 


(Continued  nn  Following  p 
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S.' 


tMa+Mb-M;-Mb  0  Ma+Mb'Ma'Mb 

x  Fx‘(jbja;LaHaLbMb;jajb)  ’ 


F  ( i  ‘  i  *  ;L  M  L.M.  ;  j  jh) 
X'ba  a  a  b  b  a  d 


(1.5-2) 


yielding 


...  is  11  b  'l  of  ki-C.  ’ ^  Q(kqq')~/2  a(W  ) 

o“(*njajbq|rnM-.,'q')k  »  — la-Tj  11(kU)  qQ 


«Vb»y;)‘-,,‘,'q  ("1>l 


k  t  i‘\ 

ooo 


x  [  a2(X) 

XM. 


+  (-0 


fq  q’  kl{(-Dq,q'  (-nt+e‘  ^  -|Wt>) 

o)1 

n 

f(Ua3b^  Ha  q'  -"a’ia'ib1} 


k  /.  .)'.+  £'-q+q’  (_-,^a+jb+ja+jb  a(qq. )  a_1(jb)  a2(j^)  a(ib) 


(-1)"  H) 


L  L.M  M.  XX’  m 
3  b  d  b  m'm" 

UiLKK 


;  i)lb*L^Ma*Nb«a*"6  o(LaLb, 


q+q' 


x  a2  ( XX'  )-i 


jJb  Lb  ■'a\Jq  La  La 

,(  Lb 

Lb 

q1' 

q*  k  q  ' 

W  j;  Ljbt  j,  ja 

k 

Mb 

Jff 

-m"  m  m* 

Q 

kal 

X  )• 

1‘ ' 

A 

X1 

'I 

k 

i  V' 

-m’ 

-M’ 

a 

Ma 

0  0 

0 

M  +M. 
a  b 

a  b 

-mj 

i”a+Mb“Ma-Mb 

0  -m 

(Equation  continued  on  following  page.) 
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drr  |  dg'  (g*  V  7^  kV .  ,(y' 


FX^J'bJa;LaMaLbnb;ja4^  F.X'  ^JL-  V’LaMaLbMb;J  aV  '  f  ‘ 


The  F,  quantities  are  expressed  in  terns  of  the  1  -1  ed 

scatterinq  matrix  usiny  (1.3-5).  This  allows  (t.b-  }  to  to* 
rewritten  upon  explicit  evaluation  of  certain  stm 


a,,(J.njajbq|i,n1j;j;q1)| 


h-  '  i  0,0  3 

(3krj: 


( k  V. >■  ‘  )  /  ►  qq  '  7  t  ■ f  '  ) 


■  1  z* t Ja jb I J ) 5 1 iM') •  K I >. •  0.10 5: 0 ■ ;  I  ,  •■(.)  ;  '  ■„  :  y 


fk  t  V 


!  i 1  •; 1  l *,  ^  ^  .  1 ,  *  -i  >  \  / ,,,,  ■  \  (  _  ;  \  I ' n  '  '  ' 


‘■‘(jj.it,)  ‘(qq')  ( - i ) 


dY  /g  ■V(Y)Rn*.'l(v)  ' 

0  0  Q  J  ay  n>.  n  ,, 


’?>  q: 


p;  qi i 

*  t  f  ( \ , 


(0  p  -aj 


|q  q'  kt 


iMa  ‘"a  ^ 


«  J  «WS<JW1iSA»Di.  (V'7l1  V 

1  a  a 


+  (-1)q+q,+  f,'+  f':  ^(j;U'iSaSb';  M  (t 


a  !V 


(Fquation  continued  on  following  pile.) 
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M  ISkgT 


;(k>.  '  yl/-  !'(kqq‘  )"^  u(Pf ' ) 


6(jajb|j^Ja)o(vqq,k|t'nOO)a'1(£)  J  «-(A)  j  d,  ^(y)^. £(y) 


1x-id+Jb‘*k  ,  ,  s ,  Mq+q'+m' 

-  MVbl'W'-"  <>(Ja3b)«(w  H-1)4 


r 

„ 

(\  v 

V  1 

A  A  ‘  £'  1 

\ 

dy 

'•  "  (  .  K  , 

.. (y)  y  «?(; 

i 

id)  n  i  n  '• 

’  A '  M 

...  m3 

o 

o 

°  . 

M  +M.  0  -M  -M. 

a  b  a  b; 

V);M 


x  (_i)v+i!  +  rVMb 

%  Jb  fl| 

fj:  j;  q 1 1 

1  a  a 

1 1  :  ( A ) 

q  q'  k 

0  v  -vi 

[0  u  -pj 

1 

v,  M.  -M  -M. 
a  b  a  b-1 

k  £  f 1 


VMb  d  -Ma-Mb 


(8’’  •' ) 


m-2 


,q’ 


|dsadSbr-S(j’j;>isasb)D^M  (Sa)D-  (Sb) 


+  (-1  ',q+q  +i  4  '  f  v+‘‘  ;:*(i 1  i 1  '  !  S  S  )  D1'  (L  )  (S  )*  ] 

1  '>  '  llbV'Vb)  UvM  1  V  i.iVV  J 

a  b 

♦  ''  (-i)k  J,‘  dt>  J<!  ^  ^:{.iaja)  '*(  ^l}Jb)  u(qq' ) 


:  -  1 


'  1 1 1  (  i)Ma,Mb-Ma-Mb  (_1}Lb+L^ 


i.  l.  m  r, 

a  b  a  b 

L  'l.'H’m: 
d  b  a  b 

v  j  i  . ' :  i ' 


ff' I  V*.  ’  )  ' 


:  .  ,)jJb  Lb  Ja!jq  La  La' 

'A'|q 


da  Lb-1  l^b  Ja  '■’a 


Lb  Lb  q'  ^ 
-Mb  Mb  Mb‘Mb 


(Fquafion  continued  en  following  page.) 
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-  C(kqq 


h 


L-Vb 


yji 


JaJb 


y 


([J](q,,«*'[J]{q,,>4?  R[q)(M  )  (-1) 

=  "  J  a  L  c  Ja 


,  ,  k+ja+jb+ja+j& 


L  L.HH,  v.V 
a  b  a  b  , 


L’L'M'M' 
a  b  a  b 


**HXV)  a=(LaLb)  1  ■;  ^ 


!  i  1  1 '  ‘  r a  L '  L  } 

•M  b  b  aH  .  a  a 


<q’  Ja  L b ^  •  -1  b  'a  Ja 


j  J 


f  «-{,  Lb  «'  If  q‘ 


’b  b 

U  «„  t-vk-^  wv\>  v;i  i";41,  -Mi  "*) 


La] 


X  x'  i*  | 


0  0  0  ]  i.M  +Mb  -m;-% 


:  1 


l«,W  0  VWV 


*  J 


^b  bb 


0  vi 


>  >  ( i 


[Ja  Jb  "a  :  jJb  Ja  Lb  ’  j  dyy,  j  dg.(g.>-  ULL 
[o  v*  -v'j(  0 


J 
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dWSidSb  s^ia>iVb: 


L„ 


-t! 


.  r  La 


*  5  (jbja>‘  iSaSb)0.-M  ‘V*  V  u  1  f'J. ' ',b;  i  i 

3 


,.(*:>*  otb  .(sn*iv  . 


(1.5-5) 


As  in  (1.3-16),  the  integration  over  the  energy  del  r.a -function 


can  be  done.  Switching  the  indices,  .r  and 


l  q  n  t 
Jt’q'n't' J 


v,  f  r 

k  | ! . ,  v ,  x'  ^ '-)  ( F> ; .  i .  ^  t{. 

n  Ja"b 


B  (^bA 


6(eqq'kK'000)  R  [0)  (e  )r[?  !  (,  . ,  )  j  dyV  R  .i  ;k 

d.j  ■'  !l  > 


-  2fi( kqq '  ;(kK.'  )*  ^  i  )q"M‘  f 


'  '  'jJb 


-aJb 
jajb 


jb  1 


x  j  j  > 

L  L.M  M.  v>V  a  D 

a  b  a  b  .,7 

L'L'M'M'  u 
a  b  a  b 


a2  (L 


lb  Jbliq 

ab  a'  f  .•  Hi- 

W  Jh  LhJ'^ 


L  *  L  I  !  L  :  L.  q 
a  a  |  1  b  b  1 

(i 


I'll’  i  -  M  M'  .M 

b  LbJ  Ja  Ja  Ja  ’  -  'b  1  b  b  n 


i  '  ! 

la  ‘a 


q  k  q  1  :  g 

i  - 

M. -M7  M’+M’-M  -M.  M  - M '  j  'M'  *'  -M'  M  f  ! f)  n  f) 

D  D  d  D  d  L>  d  d'  '  cl  a  d  d >  ' 


V 


m  +m,  -M’-m;  m’+m’-m -m.  i  !m +m. -m’ -m;  0  m‘*m*  m  -m 

ab  ab  abab'<-abab  a  ba  t 


(Equation  continue  \  on  to] lowing  pane. ) 
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''a  ^a  La  ^b  ^b  Lb]  |  i*(n+l )  Mn'tl)  i 

o  v*  -V  O  v  -u,j[‘(n+'-+-/)r(n,K'i/M 


x  (8.2)'4  |  dY  /+1  e'Y?  +  j; ]  <;v;.,,d\;ds; 


■siwiw  sVw^!5;'’  ■I"5 


. '  IV  a 

d 


'  )  ;•  '  (c.  ■  ) 

t  J, 


[y2  +  t.,  +  L  -  C.  -  i  4  ]’■'  2  L'  ^  \r'  +  !  j-  *  '  y  '  l-j  '  1  j  ) 

Ja  Jb  Ja  "b  ■’  -a  -’b  Ja  Jb 


i  (-nvtij  [Ja  La]Pb  Ji>  S  <s  <v  ^  <v]} . 

up  jO  v  -vj  (0  \.  -'.i i  a  b 

( t -5-6) 


■  ^  Q  n  * 

Further  simpl  ification  of  c"  q- n-  Jj« , p  require.,  the 
introduction  of  operators  analogous  tc  those  introduced  in 


section  1.4.  This  reductirn  of  the 


the 


)  •  h  n  i 


rol 1 ision 


integrals  is  accomplished  ir.  the  next  section. 
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1.6  Reduction  of  the 


, " ,  j. ,  '  ^  Co !  1  '  ,  ion  'nteqrals 


In  this  section,  equation  .  ,  reduced  from  a  summation 

over  nineteen  indices  and  an  into  r<.  .  i.  u  m  t  twelve  angles  to  an 

eight-fold  summation  and  a  six-told  ,  :!t  integration.  The 
reduction  closely  parallels  that  cf  1  in  section  1.4, 

and  therefore  will  be  discussed  in  •••■»;  it  k-  s  detail  than  that 


c  ,  fl  q  n  t 

of  0  Irj'.'f 


It  is  shown  in  section  1.4  th,. • 


dv  »/«•«;>’  DuVsb> 


)tv- 

1 


♦  c,  -  r  -  ;  ]l/2 

Ja  “b  -'a  Jb 


x  Ln+^(y2  +  ej'  +  cj-  ‘ 

Ja  Jb  a 


i  i 1  l  i  ■  >  i 1  1  i  i , 

v  I  ,•  a  a  a  :“h  ■;  d>  La  \  n  b  /c  i 

-  H)  n  '  L  . M  (ja}  DViM.  Sb 

vp  (0  v  -vj  (0  ..  a  b 


dSdS.dS’dS'  S(j*  j*  XlS  S,  )  '  '  IS'  S'  )D  (S' )  D  (S' 

a  b  a  b  a  b  a  b  .i  n  a  b  v  M  a  u  b 


Cah(^nt)  V  H)^fa  -  “i 

vp  [0  v 


M,  i,  C  1  hi  L 


6  b  °vM  <Sa>°-5b(Sb) 

i-  -Pj  a  b 


(I. 6-la) 
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«  I 

vu 


(-D 


v+li 


!o 


CCab(i-q,lt)  MjaJb 


X 


s  ( i '  i '  > 

tV  b 


■  r  •  c  '  \ 

>JaV 


b 


;  ,6-lb) 


Once  the  above  is  util  ; /<-d  ■ 
£ 

e val uated : 


v 


t. 


!Jb 


11.6-2) 


Making 

use  of  the 

re  1  a 

{"  L1 

La!|Ja  - 

U  L, 

b;  h 

. 

!Ja  q 

* 

V  ; 

!  1 

|o  ‘ 

-■•i  l 

and  relabeling  indices  vioi-Js 
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x  l;.  /2(y2)  dSadSbdS;dS‘[(ab(u,nt)  -S(j;jb>.|SaSb)] 


y  3  /  c  \ 


°j,  (s.^e.-vH  (sb> s  <wis;mvs;>  v,-v-a;(sb> 


(1.6-4) 

At  this  point  it  is  again  convenient  to  introduce  the 
F  b( A A  *  £ 1 ;v-v  ' )  and  Fb, (q‘ jb;6-6' )  operators  of  section  1.4. 
Now,  the  L  and  L/  summations  can  be  carried  out  as  in  (I.4-8a) 

d  D 

and  (I.4-8b).  Making  use  of  the  fact  that  these  operators  are 
hermitian  and  relabelling  indices  allows  (1,6-4)  to  be  written 


..  fz  q  n  t  ]  =  f)2v  -  J  r 

(I'q'n't'J  -  4„keT 

(LJaJb 


Y  I  l  «?(x)(p.,  :.)n(ou)_1/2  ot(jt) 


X  6(J.qq,k|£,000)R[0)(cjl)R[?)(rj,)  j  Rnl(y)  R^y) 


-  2fi(kqq')‘^  ntkU')"I/2  (i)q+q,4(Jf' 


<  c  .  i  \  !  r(n+l  1  r(n’-H) 
Tn+t+y2)r(n'+£'  + 


1/ 

(8 v'T*  l  «3(i  )  «(j-5  (p..i.)([.!](q)^[j3(q))i: 

\\*  1*  a  u  J  ,Ju  “  J- 


12X  JgJb  <Vb 

Lb 

([.J](q,)^’[j](q,))  ■,  (-1)L'+J"  *2(>A')  a2(L' L,  ) 
Jb 


ac’^Mb 

vpv'p' 


.  ,  .  c  +  Ml  fA  A'  £*H£  k  £ 

(_i)a-ni  +U+|J  -v‘v  (_i)P +  Mb 


00  00  v+v'  -v 


(Equation  continued  on  following  page.) 


X 


]dY  e”Y  ’  lnV'.'(r  )  JJ  dSadSbCcab(xx*.*; 


LA 


<  siwis.vftivciv^v 

ci  0 


dsadsb[tb'(','jb”,');i  IW'IW' 

'l 


La 


^•<s;>  °-sA'sb'  €*<v 

a  b 


1 


l. 


Next,  the  sum  over  i  is  treated: 

a 


l  a3(ja)([j](q)d’[j3(^))./ 


V, 


4  4  4 


0  a*  -a 1 J 


4  q 


a  i  a  “a 


0  p  -uj  I- 


jA+LA 

(-1)  a  a  )  , 
a  d 


(i 


making  the  introduction  of  the  f  .  (q;;)  operator:,  of  section 

a 


possible.  This  allows  the  sums  over  L1,  M‘,  i,  L,  ,  M,,  and 

dd  U  D 

to  be  evaluated: 
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,1  ^  is') 

L  M  a  aba  d 

a  a 

w 


X  “!(4>D!;M'<Sa»4(Sa>' 

a  a 


L 1 M  ‘  u 
a  a 


LbMb3  b  b  1 


(8n?)  MSa-S’)  5(Sb-S‘)  . 


<  i  .6-7) 


The  integrations  over  the  angles  S*  and  S/  can  now  be  perfonned, 

a  d 

yiel ding 


C.  q  n  t 
J'q’n'fj 


fpTT 


3A 


I 

I  l  'O'CxHp  , )  :.(0...)'V 

u;jx  a  b 


k  4;ikBT 

X  6(£qq*k /jfOOO>R^0){frj-,  )R^?)(cj^)  j  d>  y  K^y)  k'n,,.(y) 

-  2W(kqq')"^  i2(kC.fi.*  l"^  ( i  )Q+q’ +il+1L'  a(U') 


-iV 

r(n+i)  r ( r i ' + 1 j  i 

rTn+tTy;fr^'-;r+7;T| 


v-2 


(8p)~C(  )  x(j‘ )(p  ,  )  a-'(r\')  l  (-i 

i  I  -i  <  u  JaJ5 


•  *p  +  |J  1  -v-v1 


■^a^b 


vuv  p 


I  -* 


AA‘ 


A 

A' 

£'] 

l  ^ 

1 

k 

r  • 

k  q  q‘ 

0 

0 

0  j 

llo 

v+v1 

-v-v' 

-V-v'  V  v'j 

([j](q’)^,Lj'i(q,)),; 


X  ^cj^)  j  dT  v"  +  1  e'1  L^,  +  //?(y:') 


(Equation  continued  on  following  page.) 
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If  dvsbt4(iri';-"-“')4i^nt>s(j;Jbiisilsb)]D’jv(sa)D;;./(sb) 


a.  6-8) 


A  final  relabelling  of  indices  produces  an  exrression  for 
i'q'n't'jk  that  close1y  analogous  to  that  of  .qv  t' Ik 


f£  q  n  t  1  .  h2TT  ^ 


H’q'n't' 


k  4llkBT 


r 


I  ct2<A)(p.,.1);;(0Uf  ^  t(.-) 


u'Ax 


JaJb 


6(£qq,k|£,000)R[0)(ej,)R[?)(ej,)  j  dyy  ^(y)  ^(y) 


-  2fi(kqq,)"V2  fi(kt£,)'y2  ( i  )W '  a(U*) 


1 


‘A 


r(n+l )  r ( n ‘ +1  ) 


(8TT2)' 


I  a2(X\*  )a( jM (p • ,  . , )  l  (i) 

u;Jh  on 

a'  8' 


.  »u+8-o ' '  8 1 


aJb 

XX' 


X  X'  V 


0  0  0 


£  k  £' 


0  8+8'  -8-3 'j 


k  q  q ' 
-8-8*  8  3*  J 


([J](q,)©q,[J3(q,)).p  R^'^e.,) 
Jb  c  Jb 


.  £+1  -y'  ,  £' +  !A  /  \ 

dy  y  e  Ln,  • ( r  / 


dVSb  S*WISaSb>t£b('’,Jb;r''>D*V,Sb>la(,|-‘l4(Sa> 


x  Eab(XXl£1 ; - 3-3 * Kab{iqnt)S(j;j'X|SaSb)] 


(1.6-9) 


69 


As  with  equation  (1.4-15),  further  simplification  of  (1.6-9)  is 


difficult  without  restriction  of  one  or  more  of  the  indices. 
Explicit  evaluation  of  special  cases  of  the  o'  and 

!0  Q  n  f  ' 

^ t q , n , t •  |<  yield  interesting  relationships  among  the 

a  ^ ^ , 5 , f .  Some  of  these  are  examined  in  the  next 
,P  9  s  L  ^ 


section. 
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(D  q  s  t  )  (a) 

1.7  Relations  Among  the  r'|p'q's't<j 

4  Q  S  t 

Chen,  Moraal ,  and  Snider  have  examined  the  : j p • q • s 1 1 1  k 

collision  integrals  in  the  special  case  in  which  one  of  the  Bpqst 

basis  operators  is  purely  translational,  i.e.,  q'  =  t'  =  0  . 

These  exhibit  interesting  relations  due  to  the  nature  of  the 


o '  ( an  j  a  J  bq  1 51 ' n  *  J  ^  i  pO )  k  an-J  ( v.nja  Jbq  |  V  n '  j.|0)k  cross 

'  ■  \ 

sections.  The  o  S,n|'a  exhibit  similar  relations  because 
[p  U  S  U  j 

the  °'  (t‘ 0  n 1 0 ) k  and  a“(l'On'o]k  collision  integrals  are 

1  2 

themselves  closely  related  to  the  cross  sections.  In  this 
section,  the  o'  and  a"  k  collision 

integrals  are  examined  in  the  q'  =  t'  =  0  special  case  and 
relations  among  the  a|^,q  ^  scalars  analogous  to  those 

of  Chen,  Moraal,  and  Snider  are  obtained. 


If  the  restriction  is  made  that  q'  =  t'  =  0  ,  equation 
(1.6-9)  can  be  simplified  using  the  following: 


£h«>j^0)  = 

(1.7-1) 

D“'S'<V  * 

6(a'e’|00)  > 

( I . 7-2a) 

k  q  Oj 

(-1 )q_0  a_1(q)  6(q | k)  , 

(1 . 7-2b ) 

-6  B  Oj 
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([j](0)©(0)[j](0)L,  =  1  , 


(I. 7- 3a) 


p(0)/  \  } 

o 


(I. 7- 3b) 


and 


l  (P-j.  )k10)(c  •, ) Rq° ^ (e  - , )  =  6(t|0)  .  (I.7-3c) 


yayb  Jajb 


Thus  (1.6-9)  becomes 


„  (t  q  n  t 
ro  n‘0 


y  f r 
Vtt  * 


4pkgT 


l  a2(\)n((m)~^  a (£)  6(  Jtqtk | £' 000) 


dYY  Rn,(r)  Rn.t(Y) 


2-->(qqO)"^  w(qU')'V-  (-i)q"L‘fc’  a_1(q)  <t(W)  6(k|q) 


“I 


T ( n+ 1 )  i‘(n 1  +1 )  |  (  ) 

Tn+t+  /jiin'  +  r+yTJ  (8  ’ 


-2 


l  u'-(Ar)  (p.,.,) 

A)  ' 


■Vb 


(i  r 


UP 

r 


>  a  1 ' )  j  i.  q  V) 

ooo  j;o  p.  - 1-.|  J 


r  e  1  Lni  /^(y2) 


dV-’bS  <W'VV 


*  ^xvv^c^t)  ^yd^\sasb) 


(1.7-4) 


A  similar  treatment  of  equation  (1.4-15)  yields  the  important 
result: 


n'4*  q  "  1 

0  t'O  n'O 


•fi.S  ".SL  s  ■  ll-M) 


Equation  (1.7-5)  can  now  be  used  to  obtain  an  expression  Mr 
the  ojj.J  Ss.5](a)  collision  integrals,  for  q'  -  t'  -  0  , 
equation  (1.2-7)  can  be  written 


(p  q  s  tUa)  =  (_))<l+a+p  S(p' ia)u_1  (a)a(q);  (<i°b) 

°(p'0  s'Oj  'n, 


,i  /  a )  ,  ?,  q  n  t 

x  [1  +  (-1)  ]  n 1  ;psp ‘ s  1  £'0  n'^ 


:(qi?.' 


(I. 7-6 


In  detail  this  becomes 


|U)  .  ,.,ia*p  -1,A>  aioMa)?^  I  0  *  <Vl 


(p  qst|a  ,  (-Da+P  rf1  (a)  6(p’  |i 
°[p'0  S'OJ 


4'lkBT  in 


x  i(q)  .  i  m 

HnSL'n' ;psp  s 


a2(X)a((t)6(?qt!r00)  f  d.r  Rnl(v) 


rin+ii iln ' )  |  ’  (8„?)-2 


2[(i)qn+?'’  »(«‘)  f(^3^r^rT.  'TyJ 


■  / ;  v*-tJ 


x  l  a2(AA*)(pVi.)  I  0) 


x*  ?‘V?  m  fc’l 


[0  0  0  jlo 


|  dq  c‘r2  ||  dvsb  **«w'Vb> 


ia(g-n)o’6(sa)c3b<u'r;-BKab,'I,,,t)  s<W-W)f  ■ 


X 


a  result  that  is  useful  in  considering  some  special  cases  r>t  the 
collision  integrals. 

For  p  =  s  =  0  ,  the  quantity  1 ; n£‘ n ' • OOp ' f  ’  's  ff  LJ,I(J  t0  t,e 


r(q) 

C.n ft '  n '  ;00p'  s  ' 


) 2s ' +p 1 


£n,00  <‘Vn'  ,p's'  '  ,q 


h  )2s '+q 


5 ( fi.n f  *  n 1  | OOqs ' )  <S ( p ‘  |  q )  , 


( 1 . 7-9a ) 


yielding 


[0  q  0  t] 
[p 1 0  s'Oj 


: (p '  ja)  6(p‘  | q)  f;( qOq)  [1  +  (-1  )Q  J  '  V' 


,[0  q  0  0 

y  n  (q  o  s'  o! 


In  comparison,  for  p  -  1.  s  =  0  ,  it  is  found  that 


(1.7- 9b ) 


vnt'n' ; 1  Op  *  s ' 


[ ,  1 24  2s ' +q 

— j  c(?n£'n'  | OOqs  1  )  c(p':qO) 


‘  .  M2s'+q  t 


C2q+1  y 


J  t/2 


Jtnl'n'  |  COq  ( s  '  - 1 ) )  Mp'  lq-1) 


J9.L'  I  ‘  .  M,nnol  I(q) 

(?q  +  ll|  Un,lu;  qOt’n’ ;  1  Op  ’  s  ’  ’ 


(1.7-1 Od ) 


It  follows  that 


h 

23+3  __  [i 
Tq+TTT?q+2s f  -I"  3'y]  °[q+l 


q  0 
0  s' 


(q+1) 


(0  q  0  t'l'q/ 

•q  o  n, 


=  -2 


2(1-1 

V=  (  1  q  0  tl 

2q(V*l7 

L 

° ( q- 1  0  S’+1  0 

(q-D 


(1.7-12) 


analogous  to  CMS-43. 

Other  interesting  relations  are  obtained  by  t:u>v, i during 
the  case  q  =  q1  =  t*  -  0  .  Equation  (1.7-6)  becomes 


p  0  s  t 
p1  0  s'  0] 


(a 


-1 


:/ 


(-l)a  +  p  A(p':a)a"'(a)  1  (.1  +  (-1) 

v.n 
2.’  n ’ 


X  1 


(0) 


. .  1  k  On  tl 


v  n  t 1  n 1  ;psp 1  s  1  ‘  On1  0jfl 


:(0) 


Examination  cf  =  '  ,  ,  ,  ,  reveals  that 

'  o'  n  ;p‘ps  ’ 

'/ 


l(0) 

v n t. 1  n ‘  ;psp‘  s 


s+s’+p 


■>-'.//  '  0«. 


n'  ,n-»s’-s 


1  sis'  !’  (p*s  1>  jr.  )(2.  <1  )(?p+l  ) 

n  !  n '  !  ■  {  c+n+  '■  )  +n '  -*  >'-7j 


L 


‘A 


*  l 


_ (210) 


L  |s  -  n  4  '  2  L  j ' 


s  _  n  +  i+pp+ry 


p  ;  iY 


1.7-13) 


0  0  Oj 


6(>P|f p’)  , 


(1.7-14) 


which  leads  to 


jj.  0  s'  o]  =  fi(P’IP)  ^Pl*)  “''(a)  l 


n  +  /  nJi  ''^2  1 1/  \s*-s‘+p  fs is 1 ! r(p+s+  /?)  r(pfr> ' f % )  (?>'.h  jPp*  1 

L1  (1)  J  2  (/2)  hir?r+r -si’!  rr^nV-/, )7fi;r, 

r  _  (2L+1  ) _ fp  *  L|  t.,(i  0  n  t’1 

L  |s  -  n  !  r  s  -  n  +  3^~  — L  q  0  oj  I  i  0  n+5 '  "s  c  <n 


C J - 7-15) 


^  0  n_  tl  =  hVlt  ,(n 

iOn-s  0  n  4pkDT  lU) 


I  «2(X)  o  ( 1 1 0 )  j  dYY  ^(Y)  Vs.t(Y,j 


r(n+i+^,)r(n-s+l+3/^) 


f  \  i 


(8tt2  )  ~  £  ft 2  (  A  A  '  )  (  P  ■  ,  •  .  )  I 

i  '  -i  1  n  ^  h  i 


aJb  (0  0  f 


X  |  dr  /*'  e-’'2  t£V)  ||  dsa  sb  s‘(j^v:sasb) 


Eab(t.O)  Cab(£0nt)  Slj^IWi 


X 


(1.7-16) 


Equation  (1.7-15)  can  now  be  rewritten 
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a[in;st]  =  a2(i)[r(  ?>n+3£)r(i+n-s  +  3/2)]'1 

*  B 


Cl  a2(A)6(t|0)  |  dY  y2i+1  e'^  L^Vr2)  L^(y2)] 

X 


n-s 


\'  i 


a  b 
XX' 

x  (8n2)-2 


l  »2(XX')(P  )  f  d,  e'X2  l!^(v2) 

j'j'  JaJb  [O  0  OJ  i  n~s 


|  <iSadSb  s  (jjjp''  lsasb)Eab(t,o)caB(;Oiit)s(j^jjix|sasb) 


(1.7-59) 


allows  (1.7-17)  to  be  written  in  the  simple  form: 


,rfp  0  s  t 
J(p'  Os’  0 


(a) 


=  6(p*|p)6(p|a)(1/2)s+s,+p[s!s'  !r(P+st-^)r(p+s,  +  ^)]1/2 


x  [  G(ps;in)  cr[in;st]  . 
in 


(1.7-20) 


The  usefulness  of  the  above  expression  is  apparent  in  the 
treatment  of  the  following  example.  The  scalar,  3(1010)^  » 
can  be  written 


°(m2)(1)  =  %  T(7A)  l  G(ll;in)  a[in;00]  ,  (I.7-21a) 

in 


with 
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Then,  since  o[00;00j  -  0  , 

o(!o!o)(1)  =  ^  o[01;00]  +  l  o[20;00]  .  (I.7-21c) 

Similar  treatment  of  0(2000)^^  and  c:(ooio)^  yields 

o(§S8B)(2)  =  i!  oC^0;003  (1.7-22) 

and 

o(ooIE)(0)  =  |o[01;00]  ,  (1.7-23) 


from  which  it  is  apparent  that 


o( 


1 0 1 0  \  ( 1 ) 

1  L  1  0  ' 


5  -ooio- (0)  2. 

6  3 


Q  v  ^  0  0  . 


(2) 


(1.7-24) 


This  is  analogous  to  Cf‘S-4  7. 

Other  special  cases  could  be  examined  to  obtain  other  relations 

among  the  c;  |P,  ^ ;  collision  integrals.  The  examples 

l"  1  s  l  j 

treated  in  this  section  are  repr esentati ve  of  the  types  of 
manipulations  that  can  be  done  once  certain  of  the  indices  of  the 
collision  integrals  are  restricted  in  some  fashion.  The  relations 
that  can  be  found  make  it  possible  to  evaluate  collision  integrals 
based  upon  information  obtained  in  the  detailed  treatment  of  other 
collision  integrals,  at  a  great  savings  in  time  and  effort. 
Furthermore,  though  the  relations  of  this  section  apply  to 


■uaaaiilii 
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single  component  diatomic  collision  partners,  the  results  for  the 

a,(i'0  n'O  k  and  p'o  ^  collision  integrals  may  be 

generalized  to  other  collision  pairs.  This  becomes  important  in 
the  later  portions  of  this  thesis. 
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1.8  Summary 


In  summary,  the  first  section  of  this  part,  essentially  a 
review  of  Hunter's  v/ork,  provides  the  motivation  for  the 

.,.4e  Mt]  relative 

■n't'Jk  ,eldtlve 


reduction  of  the 


t  q  n  t 

rq'n't' 


and 


(e  q 


momentum  collision  integrals  of  the  remainder  of  (’art  I.  In 

section  1.2,  the  relationship  between  the  relative  momentum 

collision  integrals  and  the  kinetic  theory  cross  sections, 

a'Unj^qlf'R'j^i'q'^  and  a"  ( ?njdJbq  |  f  n‘  j  '  jfl'q'  )k  ,  of 

Hunter  and  Snider  is  established.  In  that  section,  the 

o  scalars,  introduced  in  section  1 . 1 ,  are  expressed 

p  q  s  t  J 

in  terms  of  the  relative  momentum  collision  integrals.  In 
sections  1.3  and  1.5,  the  connection  Detween  the  work  of  Hunter 
and  Snider  and  that  of  Curtiss  and  co-workers  is  laade.  In  these 
two  sections:  the  relative  momentum  collision  integrals  are 
expressed  in  terms  of  the  reduced  scattering  main*, 

b'  1  ~’a^b y  ‘ 

Sections  1.1,  1.2,  i.J,  and  1.5  set  tnr  groundwork  necessary 
for  ttie  important  result;  of  sections  l.T  and  !  6.  The  new 
resul ts  ol  Part  I  are  presented  in  these  two  sections;  the  angle 
derivative  operators,  C  ‘  aty 

E ^ ( q * j^ju1 ) ,  and  7a(q-u)  are  introduced  and  a  number  of  the 
summations  and  angle  integrations  carried  out.  In  this  manner, 
o'  f i |.  is  simplified  from  a  sixteen-fold  summation  and 

twelve-fold  angle  integration  to  a  summation  over  eight  indices 


Uqnt),  r,h( WV  ;-H'-h),  Eg(q' j^-u' ), 


q  •'  t 
■i'q'n't* 


is 


and  an  integration  over  six  angles,  while 
reduced  from  a  summation  over  nineteen  indices  and  an  integration 
over  twelve  angles  to  a  similar  eight-fold  summation  and 
six-fold  angle  integration. 

It  is  shown  in  section  !.7  that  the  rest 
q1  =  t'  =0  leads  to  a  number  of  interesting 
relative  momentum  collision  integrals  and  the 
scalars.  One  of  the  most  important  of  these 
and  a" 


tion  of 

nations  among  i.ne 
V  0  *-  t  1(a) 


P’q's't'j 

equal i ty  of 


,17  q  n  t 
[t'O  n 1 0 


■  q  n  t! 

; '  0  n  1 0  j  k 


The  simplifications  or  r  V.),],  un< 

i  v.  q  n  t  J  k 


,  '•  g  n  t  1 
I  n-  f  jk 

1  toward 


in  sections  1.4  and  1.6  provide  a  major  first  s 
calculation  of  the  transport  properties  of  a  dilute  sinqle 
component  diatomic  gas  in  the  presence  of  an  appl ied  magnetic 
field.  Though  within  the  realm  of  possibility,  such  a 
calculation  proves  extememly  difficult  due  to  the  dependence 
of  the  intermolecul ar  potential  on  three  orientation  angles  for 
a  diatomic-diatomic  collision.  On  the  other  hand,  the 


interaction  potential  for  an  atom-diatom  collision  is  a  function 

of  a  single  angle.  Calculations  of  collision  cross  sections  for 

atom-diatom  collisions  and  transport  properties  of  a  binary 

1  3 

atom-diatom  mixture  have  recently  been  completed  by  R.  Wood  for 
a  gas  of  argon  atoms  with  a  small  molecular  Mitogen  component. 
This  calculation  involves  a  gas  mixture  in  the  absence  of  a  field 
The  remainder  of  this  work  develops  the  alqehf.iir  groundwork  for 


the  calculation  of  the  transport  properties  for  a  similar  system 
in  the  presence  of  a  magnetic  field. 

In  order  to  carry  out  such  a  calculation,  it  is  first 
necessary  to  generalise  the  development  presented  in  Part  1  to  a 
binary  gas  mixture.  This  is  accomplished  in  Part  II.  Then,  in 
Part  Ill,  the  general  binary  gas  mixture  is  restricted  to  an 
atom-diatom  mixture  and  expressions  for  trie  spherical  components 
of  the  viscosity  tensor  in  the  presence  of  an  applied  magnetic 
f i el d  are  obta i ned. 


r 
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Appendix  f.A  Explicit  Expressions  tor  I(qp2;K j' u .u?‘ )  and 
/ a ( ‘li- 2 )  wi  tn  q  0,1  ,  and  2 

The  I (qpp ;tc j^nin2 '  )  sums  defined  by  equation  (.1-4-1  Oa )  are 
evaluated  using  recursion  relations  among  the  3-.;  rei  f icients ,8,9 
Their  values  for  q  =  OJ,  and  ?  are 


HOO-.Kj’VijIin'  ) 
d 


O.A-l) 


a 


■i.  .  ,  rl  >■  (•  -  )«■■)..  1  •:  1  .•  i« !  * 

M.>  •"h  :  r  1  -  ,  -|i  1  (1  1 


( T  -  A  -  2 ) 


I  ( 1 0 1K  j  '  U  i  n  •) 1  )  -  0  , 


( 1 . A- 3 ) 


Kl.-ljKj^lW) 


V.-p.-I  x>"i)  -  V.-«,  b(J>')  ■ 

(I.A-4) 


1(22;<j'yiu')  =  ^  v  +2 


*_(-  ui)  . 


(I.A-5) 


1(21  ;<  J  '  iJ  1 M 2  '  ) 

d 


"  1/26, 


( 1 .A-6a) 


'//  Kli  »*.«■)  . 


( I . A-6b) 
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II.l  Gnnrr.il 


to  H' re¬ 


used  to  obtain  • 


'<■  Ai< I  I!. 


•Y  .AS  1  XTiJKL  IN  AN 


i ■)  tv  >;,  r;r  field 


ii.v.r i/ed  Wa  1  dmann-Snider  Equation 


1  2 

ivn  1  .-i v.ann- Snider  equation  ’  is 
properties  of  a  single  component 


diatomic  gar  in  a.  ..relied  magnetic  field.  The  starting  noint  of 


the  development  > 


-f  ^  (<«+  i  L )  <t- 


(II. 1-1) 


where  ••  « ‘  r 


Bo  1  tzmann  > '  i  ’  r  !  .  * 


in  the  equilibrium  Maxwell  - 

2  3 

•nr,  I’enr-ity  operator,  normalized 


such  that 
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ana 


and  ; 
expanse : ; 


In  the  a b.. 
col  1 1 5  ion 

introo.;.;  ' 
lor?1. 
ve; or  i  tv , 
dens i . > 

•  n  th-.  a- 

int°rr  :1  ■ 
tak'1  j  tv 


r,r  ,t)  =  Tr  I  dp  f 


(0) 


( 1 1 .1 -2a ) 


c v5 ' r ,t)  -  :  r  j  dp  ,  (II. 1 -2b) 


v>'  >  >- 1  ~  o 


(II . 1 -2c) 


1  -r  I  j  ■  (  0  )  p(T1  /  ^  n  i  ill  1 

-  Tr  ,  dp  t'  [?i-  -  v:)2  +  H  ]  . 


..rsation  operator  occurring  in  the  perturbation 
or  distribution  function-density  operator, 


f  ■  f(0)  (1+4+...)  .  (11.1-3) 

v-e,  '!  and  l  are  the  linearized  Waldnann-Snider 
super. -v-erator  and  the  Larnor  precession  superoperator 
I.  T r.e  quantities  n,  Vq  ,  and  are  the 

;  v--  saonii.  number  density,  mass-average 
dens  i  tv,  respectively,  p  ’s  the  mass 
.  ;.vJ  *'  •:  s  the  internal  state  Hamiltonian 

•  •  a  .  field.  The  traces  are  taken  over 
••  iilibnum  distribution,  ,  then 

[2-nk  T]-V2e.xp[-W2.^]  ,  (II. 1-4) 

a  ,bt 


90 


where  m  is  the  molecular  mass,  W  =  (2mkgT)  ^(p-myo)  is 
the  reduced  peculiar  velocity,  and  Q  is  the  internal  state 

partition  function. 

The  general ization  of  equation  (II. 1-1)  is  accomplished  in 

direct  analogy  to  the  generalization  of  the  single  component 

4 

Boltzmann  equation  proposed  by  VJang  Chang  and  Uhlenbeck  and 

3 

discussed  by  Monchick,  Yun,  and  Mason.  In  the  absence  of  an 

external  field,  the  linearized  Waldmann-Snider  equation  for  a 

3  5 

binary  mixture  is  written  ’ 


,(0) 
'  k__ 
Tt ' 


where 


fk 


\  • k 


=  -f 


(0) 


l  '“'lA'V.1 
£  =  1 


(II. 1-5) 


t , 
k«.  k 


(2-)4tCTr,  |  dp^0’  [{ 


,dEidEk<>‘ktSkt|t|l'kkSkt>H 


»kA,h<t)t7iuM9kl>«<ek*e,-e;-e;»  (II-,-6a> 


+  (2*  i )  ^''kkEkk ‘‘kk^kk^  *k  ^k^kk^kd  *  ^kiEkkf^ 


91 


and 


k . 


-(20‘fi-Tr, 


dptf‘0)t( 


■  ^k;9kV'fE)tti“u9kl>«(Pk*erPk-p;)}  ( ii .  i  -6b) 

f  (2.7)  ','lJkf9ki  '  C  '‘‘kt^k?’^  ■  *t'llk£9kdt  luktik£>^  ‘ 


a"d  9kn 


In  the  above,  t  is  the  transition  operator,  q^1 
are  che  relative  velocities  of  particle  k  and  particle  l  before 


and  after  collision,  and  nk>,  is  the  reduced  mass  of  particles  k 


and 


Addition  of  (II. I -6a )  and  ( 1 1 . 1 -6b)  yields 


a,  +  * 

kt  k  Ki  ; 


Tr, 


dP.f  [°]U 


- 1 


dW'VSkJtK*Sk?<«'i>;> 


-<ukt3k',i«(E)ttipk,3k?*(ekterBk-e;>>  , 

*  T2T1T  <l,k(Skdt|''kt9kf>**k+fd'^kt^ 

”  <“kt9kiltd»kk9k.»]  . 


the  right  hand  side  of  which  is  clearly  the  general ization  of 
(lif  of  Ref.  7  to  multicomponent  systems.  The  essential  difference 
between  My. : k  +  and  <ii| >  lies  in  the  fact  that  in  a 

multicomponent  system,  collision  partners  may  be  of  different 

species . 


L 


The  extension  of  (11.1-5)  to  include  the  effect  of  an 
external  magnetic  field  is  immediate.  For  a  binary  mixture, 
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if'0’  K 

_J<— ,  J<  . 

>t  mk 

(n.i-8) 

g 

where  1.^  is  "he  Larmor  precession  superoperator  for  species 
k  : 

Lk*k  =  t(Jz)k’?k]  •  (H.l-9) 

In  the  following  section,  equation  (II. 1-8)  is  examined  in 

detail  and  tensor  equations  are  developed  for  the  perturbation 

operator,  ,  of  a  binary  gas  mixture.  This  allows  the  tensor 

equations  for  the  transport  properties  to  be  obtained  in  section  II. 3 

Insertion  II. 4,  seal ar equations  for  the  transport  properties  are 

obtained  from  the  tensor  equations.  These  are  written  in  terms 

of  coilisional  scalars  analogous  to  the  o(p,  ^ ,  £.1^ 

IP  *1  5  t  J 

scalars  of  Part  I.  The  new  scalars  are  expressed  in  terms  of 
binary  mixture  relative  momentum  collision  integrals  in  section  II. 5, 
drawing  to  a  close  the  formal  treatment  of  a  binary  gas  mixture 
in  an  applied  magnetic  field. 


<0)  -  f  HMC  I. 


k  L,lVk' 
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II. 2  The  Perturbation  Operator  of  a  Binary  Mixture 

The  generalisation  of  the  linearized  Waldmann- Snider 
equation  to  binary  qas  mixtures,  equation  (II. 1-8),  can  also  be 
written  in  the  compact  form 


-v, 


Uk<k  + 


2 

V 

e=l 


(<?■  4>J 

v  ki/k  k£/£ 


(II. 2-1) 


Equations  (1 1. 1-8)  and  (1 1. 2-1)  are  obtained  from  the  treatment 
of  the  following  perturbation  expansion  of  the  distribution 
function-density  operator,  fR  : 

fk  =  fk0)  H+V---)  •  (11.2-2) 

Here,  as  in  the  previous  equations,  f[0^  is  the  equilibrium 
Maxwell -Boltzmann  distribution  function-density  operator.  It 

3 

is  normalized  such  that 


where 


"Hr.o  .  Trk  f  f(°)  ^ 

(I1.2-3a) 

1/P  £  Trk  |  f*0)  Pk  dpR  =  VoCr.t)  , 

(II. 2-3b) 

P  =  l  nRmk  , 

(II. 2-3c ) 

k 
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and 

'MM  f1  [T  (sr'?°)2  +  Hi]  dfk  - 

( 1 1 . 2- 3d ) 

U^(r.t) +u(^(r,t)  =  U^(r,t)  , 

in  which  H k  is  the  internal  state  Hamiltonian  neglecting 
interactions  with  an  applied  field. 

The  requirement  that  n^  ,  v0  ,  and  be  local  values 

of  the  macroscopic  number  density,  mass-average  velocity,  and 
energy  density,  respectively,  demands  that  the  following  auxiliary 
conditions  be  satisfied: 


S- 

h- 

tl 

o 

i°‘  \  d?k  • 

(II. 2-4a ) 

°  ■  rr*  1 

,  fk0)  *k  tek-V»]  dek  - 

(II.2-4b) 

a  nd 

° =  rr*l 

r  (n\  w’k 

fk  *k  [Wk  +  l^T]  dpk  * 

(II.2-4c) 

where 

%  1/? 
-k  VkBT'  - 

mk  l/2  Pk 

!k  =  <2k‘T>  • 

(II .2-4d) 

Equations  (H.2-3a  -  II 

.2-3d)  are  satisfied  if 

fk  ^  '  (qr ^ C2nmjckgTJ  /  exp[-W^-k^j]  , 

(II. 2-5) 

■ 
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Qk  being  the  internal  state  partition  function  for  particle  k  . 
Using  the  equations  of  change,9’^0  the  differentiation  of  fj^  in 
equation  (II. 2-1)  can  be  done  to  obtain  an  expression  for 

*k  : 


/  n\  (  r>\  1  /2  C 


? 


(II. 2-6) 

2k„T  )/;>„  c  .  , 

+--t-[(|W2-1)-4~ ^(v-vo)  • 


+  nk  ?k'-k  '  cy  LV3"k 


cintT 


In  the  above,  the  bracket  notation,  [  ]^a\  indicates  a  symmetric 
traceless  tensor  of  rank  a,  cy  is  the  constant  volume  heat 
capacity  per  molecule,  c^nt  is  the  internal  heat  capacity 

3 

per  molecule,  c .  »  =  c  -7ckD  ,  and,  in  the  absence  of  external 
i nc  vc 

forces,  9  (f'^  has  not  been  allowed  to  depend  on  external 

forces) 


n,  n,  n.m, 

dk  =  v(-£)  +  [-£--i-*]7  frtp,  (II. 2-7) 

with 

V  d  =  0  .  (II. 2-8) 

k  "K 


The  perturbation,  4>k  ,  is  then  expanded  in  gradients  of 
the  macroscopic  variables  describing  the  system: 


-A  •  V  £n  T  -  Bk:  [vv0]^  +  n  V  C 


*=1 


kt 


d  -  0.7 
-l 


k!-!o 


(11.2.9) 
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in  which  ,  and  are  functions  of  the  local 

3 

velocity,  composition,  and  temperature.  Since  the  gradients 
in  equations  (11.2-6)  and  (II. 2-9)  are  independent,  the  integral 
equations  separate: 


2(«kwk-iH-’y)  •  UA.  I  («t,8k»«-kJt)  • 


(II. 2-10) 


where  u  is  the  second  rank  unit  tensor. 


2knl  '/, 


<Ti7>  C(Mj-f)*-V-*]Hk  I,  • 

k  B  t=l  (II. 2-11) 


c.  .  o  h'-<h;>. 

i nt  r/Mp  k  k  k-i 


C  )  "  r - r 

Cv  J  K  Cint 


’  ukV  <"ktVW 


(II. 2-12) 


2kRT  1/2  I 

n^-k^kj~5ki^  =  lLk^-ki"-kj) 


+  y  fai.'  (C.  .-C,  . )  +  47.1  (C„.-C  .)]  , 
^  L  ki  -ki  -kj  '  kii v  -  2.  l  -ij/J  ’ 


(II. 2-13) 


The  last  equation  is  obtained  by  consideration  of  the  condition 
on  the  d^  stated  in  (II. 2-8). 


Equation  (II. 2-10)  requires  that  §k  be  a  symmetric, 

traceless  tensor.  Thus,  it  automatical ly  satisfies  all  three 

8 

of  the  auxiliary  conditions  expressed  in  equations  (II. 2-4). 
The  Ak  and  automatically  satisfy  (II. 2-4a )  and  (II.2-4c), 

however,  (II. 2-4b)  yields  the  two  auxiliary  conditions: 

0  =  l  ^\Trk  J  fk°}  dPk  til. 2-14) 

ana 

o  -  I-kM  fk0)  <^k1-CkJ]-Uk)  dpk  (11.2-15) 
k  K 

The  Dk  on  the  other  hand,  satisfy  (II. 2-4b)  automatically, 
while  (II. 2-4a )  and  (II.2-4c)  provide  the  remaining  two  auxiliary 
conditions: 

0  =  Trk  |  ^°}  Dk  dpk  (II  .2-1 6a) 

and 

Hf 

0  "  jjTlk  I  fk0)  “k  ['V Sjrl  %  ■  <n-2-,6b> 

The  integral  equations,  (II. 2-10)  through  (II. 2-13),  along 
with  the  auxiliary  conditions  stated  in  equations  (II. 2-14) 
through  (II. 2-16),  may  now  be  solved  for  the  Ak  ,  §k  ,  Cki  , 
and  Dk  .  It  is  useful  at  this  point,  however,  to  obtain  tensor 
equations  for  the  transport  coefficients  in  terms  of  these 
quantities.  This  is  done  in  the  following  section. 
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II. 3  Transport  Properties  of  a  Binary  Mixture 

In  this  section,  the  perturbation  operator  of  section  II. 2 

is  used  to  obtain  expressions  for  the  pressure  tensor,  P  ,  the 

diffusion  velocity  <V|>AV  »  the  energy  flux,  g  ,  and  the 

anisotropic,  (symmetric,  traceless)  part  of  the  dielectric 

tensor,  [c]v  ‘  ,  in  terms  of  the  ,  0k  »  ,  and  Dk  . 

Comparison  with  phenomenological  expressions  for  P  ,  <Vj >AV  , 

(2) 

q  ,  and  [e]v  in  terms  of  the  experimentally  ineasureable 
transport  coefficients  yields  expressions  for  the  transport 
properties  in  terms  of  the  A^  ,  §k  ,  ,  and  Dk  . 

The  pressure  tensor  is  given  by  the  following  phenomenological 
relation:  3 


P  =  pU  -  2n  :  [WoT  '  -  <(V-v0)  , 


(II. 3-1) 


in  which  p  is  the  scalar  pressure,  [Vv0]^  is  the  rate  of 
shear  tensor  (symmetric  and  traceless),  n  is  the  shear 
viscosity  tensor  (of  rank,  4),  and  <  is  the  bulk  viscosity 
tensor  (of  rank,  2).  In  terms  of  the  distribution  function- 

9 

density  operator,  f.  , 


2  ■  2kBT  jj  Tkk  {  fk  "k  -k  d?( 


'  2kBT 


l  Trk  f  fk0)  0 V  “k  “k  dek 


(II. 3-2) 


i 

1 
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Neglecting  the  terms  involving  temperature  and  concentration 
gradients. 


E  *  2kBT£(Trk 


fk01“kVPk  -  Trk  |  f[0,“k"kik  ;  «oOpk 

(II. 3-3) 


-  Tr, 


f[0)«kW?°d<V 


In  the  above,  and  in  the  equations  to  follow,  the  integration 
over  pk  can  equally  well  be  taken  to  be  an  integration  over 
particle  momentum  or  peculiar  momentum. 

9 

The  following  identifications  can  now  be  made: 


I  Trk  F(Wk)WkWkdpk  =  Trk  F(Wk)W2dpk  ,  (II. 3-4) 


2kBT 


(0). 


ielding  p  =  -y- £  Trk  j  Wkd?k  =  nkBT  ;  ( II .3-5) 


v  l  Trk  |  <0)ykykikdPk  * 


(II. 3-6) 


since  Bk  is  symmetric  and  traceless;  and 


8  *  2k6T  1  Trk 


fk°^-k-kDkdpk 


(II. 3-7) 
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Similarly,  a  phenomenological  relation  involving  the 

3 

diffusion  velocity  is  given  by  the  expression: 

n-m  T 

Vk^Av  ■  --  l  T  ■  (''•»-«) 

K 

in  which  ,  and  0^  are  the  multicomponent  diffusion  and 

thermal  diffusion  tensors,  respectively.  In  terms  of  the 

3  9 

distribution  function-density  operator,  ’ 


<V 


k"AV 


--  Tr 
k  k 


(11.3-9) 


and  consequently,  neglecting  terms  involving  vv0  and  V-v0  , 


Vffi*  =  'mkTrk 


fk0)- A  •  J  T  dBk 


*  V  l  Trk  I  fk0)-k?kt'-i  d?k  • 


(II. 3-10) 


since 


A  I  fl0>h  A  =  0 


(II. 3-11) 


In  this  case,  the  following  identifications  can  be  made 
immediately:  ^  ^ 


-k  i 


(II. 3-12) 


- - 
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and 


»k  ■  mkTrk 


f  fk0>-k-k  dfk 


(II. 3-13) 


3  9 

The  phenomenological  expression  for  the  energy  flux  ’  is 


=  k, 


>T  [  '?♦ 


5  .  ‘^Pk 


~kjjT  "'k'-k^AV 


(II. 3-14) 


nkcT 


1  J 


Dl  *  d.  , 


6  t  w y k 


in  which  a3  is  the  thermal  conductivity  tensor  for  a  gas  mixture 

of  uniform  composition,  that  is,  \0  is  the  thermal  conductivity 

tensor  in  the  limit  that  no  thermal  diffusion  has  occurred. 

If  equation  (II. 3-8)  is  solved  for  the  d^  in  terms  of 

<Vk>AV  and  VT  »  and  the  result  inserted  in  (II. 3-14),  the 

new  coefficient  of  VT  is  -X^  ,  the  thermal  conductivity 

tensor  which  is  ordinarily  measured. 

In  terms  of  f^  ,  the  enerqy  flux,  neglecting  terms 

3  8  9 

involving  Vv0  and  V-v0  ,  is  given  by:  ’  ’ 
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„  2k  T  2 

3  '  i  ‘if1  rrK 

k  k 


fk*kB™k  *  Hk'-k  dSk 


I(^1/2<-Trk|fkQ><kB™k*HP“k 


+  n  7  Tr, 


l  fw 


k  '(VWk  +  Hi>“A;.-3k 


since 


f  f(0) 


"i>“k  (?k  *  0 


Now,  using  (II. 3-10) , 


c  <«k>k 

k8T  l  (?+_kgT)r1k<^k>AV 


,  <H!> 


X  (Tr, 


'f(0)ykAk.v«nTdPk-2nTrk[  f(0,VkCM-d 

?  J 


allowing  (1 1. 3-1 5)  to  be  rewritten: 


r  <H!>. 

3  *  V  {  (7*”k^)"k<!!k>flV 


2k  T  l/2  ,  k  H!  ,  <H!> 

-kBT  (-i T> 


Hf  <H' 

I  Trfc  f 


Ak  •  7  InT  dpk 

dpk}  , 

(II. 3-15) 

(II. 3-16) 


(II. 3-17) 

t  V  * 

(II. 3-18) 

rr^k*  •  V?  k} 


By  using  equations  (II. 2-11)  and  (11.2-13)  along  with  the  symmetry 
6  8  12 

properties  ’  ’  of  ,  and  L k  ,  it  is  possible  to 

show  that  (see  Appendix  1 1. A) 


2kRT  */2  f  (01  Hk  5  <HL>k 

l  (  — 5_)  y  Tr.  fJ0,(Wt  +  A-|--tV)wlCu.  •  d0  dp 

£  '  mfc  £  k  J  k  '  k  kgT  2  k^T  -k-kt  -i  - 


-  >'  nV  Pi'?,- 


(II. 3-19) 


£'  £ 


yielding 


,  <h:>. 

kBT  £  h  +  Yg'HV-k^' AV 
2kn7 


k  "k 


wo  (nv  H’  r 

kBT  l  (ir)  Trk  T)^  +  CT-f--xr^k*^nTdPk 


<h:>, 


nkRT  T  DT  •  d  , 

B  “  m  n  =k,  -  £ 

£  f  £ 


(II. 3-20) 


from  which  the  identification  can  be  made  that 


2knT  j /;  /n\  c 

kB  i  <-5®->  Trk  | f  k  1  <Mk  - 1*  -  A  d?k ' 


(II. 3-21) 

Finally,  if  only  gradients  in  v0  need  be  considered. 


the  symmetric,  traceless  part  of  the  dielectric  tensor  is  given 

13 

by  the  phenomenological  expression: 
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OJ(2)  =  -2|  :  [7v0J(2)  ,  (II. 3-22) 

in  which  §  is  the  flow  birefringence  tensor  (of  rank,  4). 

1 3 

In  terms  of  the  distribution  function-density  operator, 


[|](2)=  -2:r  &  I  nk  (ot„-  ct,  )kTrk  [  dpk  fR  [j'UJ-J)}’ ^^k^ 

~  k  ii  L  ■  . 

*  2*  I  "|it»||-’i>kTrkj'  05kf(0,lj((4-I»"/2UkJkl(2|>lk:!Yo  • 

(II. 3-23) 


in  which  ay  and  are  the  electric  polarizabilities  of  a 

molecule  parallel  and  perpendicular  to  the  symmetry  axis  of  the 

molecule,  respectively,  and  in  which  the  second  rank  tensor, 

{/^j£(jj^)]"1/2[JkJk]^}  ,  characterizes  the  angular 

14 

momentum  polarization  of  the  molecules  of  species  k  . 
Again,  the  identification  is  immediate: 


i  ■  "k(»in>k  Trk  |  §k  d2k  • 

(II. 3-24) 
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Using  the  perturbation  of  section  II. 2,  expressions  have 
been  obtained  for  the  shear  viscosity  tensor,  n  ,  the  bulk 
viscosity  tensor,  <  ,  the  multicomponent  diffusion  tensor, 

,  the  thermal  diffusion  tensor,  ,  the  uniform  composition 
thermal  conductivity  tensor,  ,  and  the  flow  birefringence 
tensor,  §  .  In  the  following  section,  scalar  equations  are 
obtained  for  the  above  quantities,  yielding  expressions  more 
useful  for  computational  purposes. 


■  J£* * 


II. 4  Scalar  Equations  for  the  Transport  Prooert'es 


Equations  (II. 2-10)  throe:-:  (II. 2-13)  can  to  written  in 
terms  of  tensors  in  the  Wang  Chang-Uhlenbeck  basis.  Expansion 
of  the  unknown  tensors  of  the  perturoation  operator  appearing 
in  these  equations  over  a  total  polarization  basis  leads  to 

g 

a  partial  uncoupling  aue  to  the  conservation  of  oolarization. 
This  yields  the  desired  scalar  equations  for  the  transport 
properties. 

In  Ref.  12,  the  Wang  Chanc-Uhl enbeck  basis  elements  are 
defined : 


?k,pqSt  =  LF5(«,;)[,'i!](qlR<,,)(Hk/k3T)  .  (II.4 


the  {.^(W^)  being  velocity  tensors  composed  of  a  product  of 
an  irreducible  Cartesian  tensor,  [W]^  ,  of  weight  p  ,  and 
an  appropriately  normalized  ^  associated  Laguerre  function, 

L Cp+1A ) (IJ2 )  _  The  left  hand  sides  of  equations  (II. 2-10) 
through  (II. 2-13)  are  expressible  in  terms  of  the 
-k  pqst  e^en1en,s  Just  defined.  (II. 2-10)  through  (II. 2-13) 
can  then  be  written 


iLk!k+  J,  =  2(-k-k '  rk=) 


»  /2  L20(Wk) 


(II. 4-2 


=  /2  B 

-K.2000 


2  2k  T  1/?  5  HI 


L  4S'"R'  a*  kB  '  k  b 


/~5  n  +  )  i  nt  p  -i 

'  mk  *-  J  1  =k;loio  \  kg  =k;iooi^  ’ 


(II. 4-3 


ukDk + mv 


t-/|  Rl0)  (Hk/ k  T)] 

cv  cint  «  kB  k  b 

Cint  r  /Tr  kB  C  int  p 

c  L~  3  -k;ooio~  \  kD  =k;000iJ  ’ 
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uk‘Skr?kj>  *  j,  «‘lkt‘SM-Skj>,St‘£.i-S«>’ 


2kRT  i j2  i 

<-S7>  K^k<W 

2k„T  1/2  l 

'  7r^fc10(“k>(v5ki: 


(II. 4-5) 


kqT  1/2  I 

§k;iooo^5kj-<5ki  *  ’ 


The  transformation  of  (1 1.4-2),  (1 1-4-3),  and  (1 1.4-5)  to 
scalar  equations  is  accomplished  by  writing  them  in  terms  of 
their  spherical  components.  As  in  Hunter's  treatment  of  a 

g 

single  component  system,  the  spherical  basis  chosen  is  that 

of  Chen,  Moraal ,  and  Snider.  12  Their  basis  tensors,  epm  , 
transform  under  rotation  according  to  the  irreducible 
representation  of  the  rotation  group.  Furthermore, 


6 


m,m' 


(II. 4-6) 


where 


and 


&  -  ‘fK>* 


m 


Pjw  = 


(II. 4-7) 


-  V 
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where  is  the  projection  operator  onto  the  space  of  p-th 

rank  symmetric  traceless  tensors.  Consequently,  the  are 

orthonormal  and  complete. 

The  Bk  appearing  on  the  LHS  of  (II. 4-2)  are  symmetric 
traceless  tensors  of  rank  2,  and  consequently,  their  spherical 
components  are  obtained  by  dotting  them  with  e2m  : 

Bk  =  e2m  ©2  Bk  .  (II. 4-8) 

The  and  C^.  of  equations  (II. 4-3)  and  (II.  4-5)  are  of 
rank  1,  and  so  their  spherical  components  are  obtained  by 
dotting  them  with  e^m  : 


and 


(II.4-9a) 


0‘ki 


( 1 1 .4~9b) 


Note  that  (11.4-8)  and  (II . 4 - 9a )  are  simply  component  versions 
of  H-44  and  H-52. 

Using  the  above  and  the  fact  that  the  spherical  components 
of  the  §k;post  al"e  9iven  by  ° 


e 


mo 

k;post 


J5"1 


©K  B 


-k;post 


(II. 4-10) 


II  03 


no 

it  Is  possible  to  obtain  scalar  equations  from  (II. 4-2), 

(II. 4-3),  and  (II. 4-5)  by  dotting  (II. 4-2)  with  e^m  ,  and 
(II. 4-3)  and  (II. 4-5)  with  e'm  .  In  addition,  noting  that 
is  a  zeroth  rank  tensor,  i.e.. 


B  =  B 

=k;oost  k;oost  ’ 


(II. 4-11) 


(II. 4-2)  through  (II. 4-5)  can  be  written  in  scalar  form: 


/2  Bm° 
c  k;?oog 


u/k 1  l  I'WW  • 

i  =  I 


(II. 4-13) 


,mo 


'int 


/kck 
*B  int 


(II. 4-14) 


and 


,kB^\*/2  1  Dmo  /x  ,  »  _  . ,  ,.m  rm  . 

mk  nk  Bk’  1000  '  k-j  ^ki  ~  ^k  ^ki’Ckj 


(II. 4-15) 


*  J,  ■ 


Ill 


The  spherical  components  of  the  Wang  Chang-Uhlenbeck 
basis  operators  are  given  by 


k;pqst 

=  e1*  ©p  B,  t  0q  eqv  , 

=  =k;pqst  = 

( I 1 .4-1 6a) 

which  in  turn  yield  the  spherical  components  of  the 

total 

Q 

polarization  basis: 

D  (  3  )°t 

Dk;pqst 

(-i )a+a  v  |P<|  a)  b[;v  ,  . 
po  v  -ctJ  kJPqst 

(II. 4- 1 6b) 

The  components  of  the 

perturbation  operatcr  are  now 

expanded 

over  the  total  polarization  basis: 

■ 

r  BPqst  (a)m 

pqsta  (a)m;t  *iPqst  * 

(II. 4-1 7a ) 

■ 

y  Apqst  B(a)m 
pqsta  (a)m;£  £;pqst  ’ 

( 1 1 . 4-1 7b) 

\  - 

and 

n0  .  v  npqst  R(a)0 
'  "pqsta  W*  • 

( 1 1 . 4 - 1 7c ) 

rm 

Si 

y  cpqst  b^3^1 

pqsta  {a)mUi  *;pqst  * 

( 1 1 . 4- 1 7d ) 

in  which  the  B„a^V 
£;pqst 

are  the  basis  elements,  and 

Rpqst 

A(aRt  •  D(a)0,t  •  an<1  C5U;„  are  the  unknown  «Pai>^<>" 


coefficients. 
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Utilizing  the  above, 


rc  Dk;2000 


Jl  pqsta  {[BU)m;k(^k2+1:kiLk)Bk;pqst] 


4-  fR^S^  1) 

( a  )nr, 2  KZ°&;pqst  ’ 


(II. 4-18) 


•  „k 


,kBT|'/-r  Cirn  .(l)m  ,-5.(l).  i 

'  m.  \  k„  k ; 1  -  - 1  ^  ?  k;!on 


(II. 4-19) 


1=1  pqsta  v 


cint  r  /T„  <^3Cint  „  .  .  9fU 

c  ^  J8k;00l0  c.  Bk ;  0  0  C  i  ^  (II. 4-20) 

v  int 

•  !,  l  +  »B)0;  >  ' 

1=1  pqsta 


and 


fVV7*-  _L  B(1)m  u  -6  ) 

'  Bk;ioao  5kj  5ki ' 


J,  Jta  ««$!&*«  -  cU)™;y)(Vi5Ulk>Bktpqsti  <"-4-2’> 


,(a)m 


H 
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The  basis  elements,  B[a^ac. 

kipqst 


are  normalized  such  that 


«r<rR  (a)a  |  a  (a '  )<*' 

<<B  k;pqst* “kjp'q' s't‘  >>: 


6(aa|a'a' )6(pqst|p,q's't' )  , 

(II. 4-22) 


where  the  inner  product  for  mixtures  is  defined 


12 


as 


5  Trk  (  dBk  f(0)  Ak3k  •  (n-4-«) 

Multiplying  equations  (II. 4-18)  through  (II. 4-21)  from  the  left 
by  and  taking  the  inner  product,  the  following 

set  of  matrix  equations  is  obtained: 


ft  ^p'q's't'a*  [20002) 


?  r  ;rBPqst  ^ „(a ’  )m'  L,)iB^a^  ».] 

ti,  pqka  [  (a)m;k  kSP'q's'f  kl  kt  k  k,pqst  k 

(1 1. 4-24) 


+  [B?qf  «B)“  » ,  ] : 

L  (a)m;f  k ; p  q  s  t  W  £;pqst  kJ 


.(a '  )m' 


(a)m 


koT  »/.'  jc;n,  /r 

(_5_)  [  f.fp'o's’t'a' 110011)  -  /%  ^p'q's't'a'noiOl)] 

l  Xg  *- 

£=1  pqsta 

*  raPqst  ^.^p(a’  }ir‘  p(a)m  v  r 

t  (a)m;:_  k’.D'o's’t’1  kk  1  k;pqst'  k  ’ 


(II. 4-25) 
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r~T 


— ^  [/j  5(p‘ q*  s'  t 'a '  1 00100)  +  ^ (p'q's*  t'a' |00010)] 


V  V  .[npqst  --c  r(^  ‘  )m  ’ 
"  ~  l.  I  l-UfA\n-k 


?,=  1  pasta 


U  )0 ;  k  p i q  ■  S  •  t  •  I  < - 6 ^*-1, >  I  pqs t  >>  k J 


+  ro?q^  «B^a  i if  »  ii 

U(a)0;i  k ; p 1  q 1  s ’  t *  ‘  kS.1  ;.;pqst  ’ 


(II. 4-26) 


(^'"iflp’qVtVHOOODUy.S^ 


f  Tfrpqs,t  -  rpqct  l 
‘-'L(a)m;ki  L(,a)m;kj 


t=l  pqstu 


*  <<Bktp!q:s'fl<V3Mtk)i8ktMStS>kl 


(II. 4-27) 


t[(c?s;ii-cnM;u)<<Bk-pvs't'i<,'ai8S!p;t>>ki) 


The  fact  that  and  a!j^  are  invariant  under  all 

g 

rotations,  and  consequently  conserve  both  total  polarization 
indices,  allows  their  matrix  elements  to  be  written 


«B[a'K  i  * .  lK,o  i  b[3  .».  =  6  ( a  ’  cx '  jau)«  B^?  ,  |B^a 

k;  p‘  q  *  s '  t ' 1  k«. 1  k;pqst  k  v  1  '  k;p'q's't"  kr  k;pq 


pqst  k  ’ 


(II . 4-28a ) 
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and 


(a'  )a 
k;p*q‘ s‘ t* 


(a)a  > 
k;pqst 


k 


6(a'a' | 


aa)<B 


(d)n 

k ;  p  ‘  q 1  s '  t 


■Kd6 


(a)a  . 
k;pqst‘ 


k 


(II. 4-28b) 


It  is  convenient  to  define  two  scalar  quantities  as  the  matrix 
elements  appearing  on  the  RHS  of  ( 1 1 . 4 - ?8a )  and  (II. 4- 28b) : 


.  p‘  q’ 

IP  q 


t* 

t 


(a) 

U 


V  q'  S'  f)(a) 
,p  q  s  t  j  k; 


^  n  ( a  )a  i  j  g(a  )o 

K ; p ’  q 1  s 1 1 1  1  let1  k;pqst 


» 


k  ’ 

(II. 4- 29a ) 


R(ah 

k ;  p'  q 1  s 1 1 1 


'V» 


(an 
f. ;  pqst 


> 


J.  ‘ 


(II.4-29b) 


,  on  the  other  hand,  is  invariant  only  to  rotations  about  the 

g 

field  direction.  Consequently,  it  conserves  only  the  u 

index,  and  explicitly. 


(a '  )a 
k;p’q's't’ 


(a)u 

k;pqst 


where 


(u>L)k  SCp'n's't’a'  ipqstfK)Lp“ha'a)  , 

( 1 1 .4-30a ) 


(-1 )a+a'  + ‘(-1 )P+q+1/(2a+l )  (2a ‘  +  1 ).'q(q+l )|2q+l I 


fl  a  a'j/  1  a  a'. 

{?  a  -oj  [  p  q  qj  * 

(II. 4-30b  ) 


116 


Using  the  above,  equations  (II. 4-24}  through  (II. 4-27)  can 


be  rewritten: 


2 

5 (pqsta 120002}  -  ~  ,  B??!5' t*r5  ^‘fp  q  s  t|(a)l 

;  =  ]  p'q's't'a'  ^m;k  1  aa'~  Ip'  q'  s‘  t'm  J 

*  B(a')";r[5aa's''^  g'  l  ^kil4iV\Vl^!MsOl‘;!Wm. 


(II. 4-31) 


.v4«/»  R„t 


(--)  [  Mpqstaj  10011)  -  /  |  5 (pqsta i  101 01 )] 

2 

=  y  "  aPV  S‘  t'fr  c; i  f  p  q  S  t  1(a)-, 

til  p'q's't'a'  'a‘:ir'  k  6aa'~  lp‘  O'  S'  V;kl  } 

+  A(a*)ni;5  ^aa,S"  p  q'  s'  r j  fi6k£K)k'(^lPqSt)Lpq)(aa,3}' 


(II. 4-32) 


AJ: 


3  -S  ( pqsta  1  00100)  +  — 6 ( pqsta |  C001 0)  ] 

int 


=  y  )  rDP*q' s’ t’r^  C,  p  q  S  t  )  (a). 

Jl  p'q's't’a'’  (a')0;k  aa’  IP’  q'  S'  t  *j  kJt  1 

+  ^ROa^aa'^'tp'  5'  s'  t]k,  +i6k.t(uL)k<S(r,f,s,Li,pqst)Lp°q  (aa^}* 

(II. 4-33) 


and 
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(^)’/!i6(pqsta!10001)(5kj-5k() 


r  y  {/rP'q's't'  _  rP'q's't'  c.fp  q  S  t] 

£=1  jfc^a'  L(a')miki  L(a'  )m;kjH6aa'b  [p*  q'  s'  tj 


(II. 4-34) 


The  above  set  of  matrix  equations  can  be  solved  to  obtain  the 
expansion  coefficients,  , 

and  ,  subject  to  the  following  auxiliary  conditions:® 


n.  *n.  V"'0? 
£  k  k  ln„k 


(II. 4-35) 


0  =  l  \  ^  (Clm?ki  ”  Clm;  kj  ^  ’  <II4-36> 


D?ojS;k 


(II. 4-37) 


0  =  l  r'k  '  'I  C  °(0^;k  “  ^  T  n?o]6;k}  * 

(II .4-38) 

Recalling  equation  (II. 3-6),  the  shear  viscosity  tensor 
can  be  written  (since  Bk  is  symmetric  and  traceless) 


1 
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1  „  „  (21 
n  =  kni  >  n,.  -CB' 


.'2  *B'  i  "k  '^ooo'^k  V  (H.4-39a) 


Expanding  in  terms  of  spherical  components  and  then  over  the  total 
polarization  basis  yields 


n  = 


1 


/2  B 


kJ  l  n 


k 

mm 


e*  ‘  B?'  e2  ,> 

=m  k ;  2  0  00’  k  -  m 


=  /2  kBT  l  nk  ^ 

k  mm  pqsta 


2m  D(2)iri  |  D(a)m'  ^Dpqst 


e  "  «  B 


k;:  ooo  I  Bk  ;pqst~  ”(a)m'  ;k  Sm' 
( 1 1 . 4-39b) 


Using  the  above,  the  viscosity  can  be  expressed  as 


8 


with 


n  = 


m 


l  I 

m 

kBT 

/2 


2m 


n  e" 
m  =m 


n  r2000 
k  (2)m;k 


(II. 4-40) 


(II. 4-41) 


Similarly,  since  Dk  is  a  scalar,  equation  (II. 3-7)  for 
the  bulk  viscosity  tensor  can  be  written 


r 


/2  k  T  V  n  <*78 _  1^5.  B(°)  U  )  I 
d  V  £  k  ^=k;2000  3  Bk;Ooio=  }  |Dk  ^k 


,(0). 


/2  k„T  /  n.{-5<e[°J n.JD^>.  U  +  l  e2m<B[2),T1  |  dJ°*> 

B  E  k  3  k ;  oo  l  o  1  k  k=  “  =  k ;  2  o  o  o 1  k  k 

In  ITl 

/o  k  T  7  n  fe2°D200°  -  ^  u  ) 

/2  kgl  l  ryg  D^2j0.k  -j- U  D(0j0.k}  . 


(II. 4-42) 
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Consequently,  <  can  be  expressed  as 


g  =  e20?  -  u  nv  , 


c  ^  kBT  £  nk  D?2?0;k 


r^_  l  t  v  _  pvOOIO 

*  3  V  £  nk  U(0)0;k 


(II. 4-43) 


(II. 4-44) 


(II. 4-45) 


Recalling  equation  (II. 3-21),  the  uniform  composition 
thermal  conductivity  tensor  can  be  written 


kB  l  »k,iic..»li|e>k  • 


(11.4-46) 


In  terms  of  spherical  components  this  becomes  (after  expansion 
over  the  total  polarization  basis) 


i  -  v  v  /  B  r  r  lm 

=0  ~  ^  m  nk  <-  4  - 

k  mk  mm'  pqsta 


'  2  k;1010  k_  “k:lom  I  kiDast^k  iaW:k  -m' 


kR  k ; i oo i  1  k ;pqst  k  (a)m';k  -m' 


(I I. 4-47) 


HP-1  uujuii.  i." 
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l  ■ 

L 


This  leads  directly  to  the  expression 


-  v 

£0  -  L  ?  *o  ?  , 


(II. 4-48) 


k  T  n  c 

x*?  =  ic  7  r  fl1001  J~5  A1010  1  (II  4-49) 

A°  “  £  "n  kB  A(l)tn;k  -  /  2  A(l)m;kJ-U1  ^ 


The  flow  birefringence  tensor  is  treated  in  direct  analogy 
to  n  •  The  spherical  components  of  &  are  given  by 


*(««»■  (n-4'50> 


The  remaining  two  transport  tensors,  the  multicomponent  diffusion 

tensor,  ,  and  the  thermal  diffusion  tensor,  ,  are  treated 

in  direct  analogy  to  A0  •  The  resulting  expressions  for  the 

spherical  components  of  D  and  dJ  are 

“  k  =  k 


nm  =  p  nk  / kBT \ l/2riooo 
ukt  nmi  '  mk  ;  L(1  )m;k2 


(II. 4-51) 


Dk;m  =  HkBT)1/2nkA(l)m;k  *  (II. 4-52) 
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The  above  expressions  for  the  spherical  components  of  the 
transport  tensors  involve  specific  expansion  coefficients, 
which  in  turn  can  be  expressed  in  terms  of  the  S'  {jj,  * 

S"  ^  ^  ,  an£*  Lpq^aa'^  elements,  once  the 

set  of  equations  (II. 4-31)  through  (II. 4-34)  has  been  solved. 

In  the  following  section.  S'  r ,  tl)a  and 

S"  ^  are  examined  and  related  to  the  collision 

integrals  of  Part  I.  To  the  extent  that  these  can  be  determined 
computationally,  the  transport  properties  of  a  binary  gas  mixture 
in  the  presence  of  an  applied  field  can  now  be  determined  and 
compared  with  results  obtained  experimentally. 


II. 5  Expansion  of  the  S*  (jj,  jj.  *,  and  S"  jjj,  jj.  * 

in  Terms  of  the  Re’  '^e  Momentum  Collision  Integrals 


As  in  Part  I,  it  is  possible  to  expand  the  scalars, 

S’fp'  q'  s’  t‘  and  S"  o'  q1  s'  t'  ’  aPPearinq  in  the 
matrix  equations  for  the  expansion  coefficients,  in  terms  of  the 


...  .  .  ,  ,  q  n  t)  v  1 9  q  n  t]  v 

collision  integrals,  o  je,  jj,  ^  and  q,  t,J  k  • 

Here,  to  avoid  confusion  with  the  k  and  l  indices  that  appear 

in  the  collision  integrals,  the  species  of  the  collision  pair 

are  labelled  by  the  indices  5  and  v  .  The  development  of 

this  section  relies  heavily  upon  analogy  to  Part  1,  and 

consequently,  avoids  details  already  presented  there. 

If,  as  in  Part  I,  only  a  single  component  gas  is  being 
considered,® 


«Ofl|BpVs;t;*  *  "<8&>,/2°(p’  n'  r  f'](al«(a«|av>  • 


(II. 5-1) 


Recalling  that  can  be  separated  into  tfT'  and  fl"  ,  this 
could  be  rewritten  as 
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in  which 


(p  q  s  t 

'  (p'  q'  s'  t’ 


(a)  _ 


(II. 5-3) 


and 


l  l  M<a>  ■ 

(II. 5-4) 


From  the  development  of  section  1.2  it  is  also  possible  to  write 
the  above  as 


,  fp  q  s  t  |  (a) 

0  lp'  q'  S'  t'j 

■  f-u^'v  (-l)k«2<k)c(kq'q)>/4;,  2'k)  l  c(k«')l/! 
r  vH  r  /  on 


i 


(k) 


,  ft  q  n  t 


tnt'n' ; psp '  s'  (>'  q'  n'  t’ 


(II. 5-5) 


and 


p  q  s  t  )(a) 
P’  q*  s’  t'j 


tn 


•  (-l)q+a+P  ^  (-l)ko-(k),(kq'q)'/J(j,  J'  *}  l  Q(kU')>/*(-l)1' 

(1 1.5-6) 


T(k)  J*  q  n  t  . 

Hns-'n'  ;psp's'  q'  n'  t'j 


k  ‘ 


Equations  (II. 5-3)  and  (II. 5-4)  can  be  generalized  for 


binary  mixtures  as  follows: 


^-0(3)®  |/ji  |n(a‘ )<*'  ^ 

^  a ;  pqst  ‘  sv '  6 ; 6 


and 


"Pr  -1/2 

ny(airf)  Sfaula'c-Jo' 
B 


p  q  s 
p‘  q'  s' 


t  ](a) 

t'Jfiv 


(II. 5-7) 


I aj"  |ola  )a  >. 
B6;pqst'  <Sv‘  & 


(a*  )a ' 


TTp  -l  /  2 

nJw^)  <s(ja,aS,)o" 


pqst 
p'  q'  S'  t' 


(a) 

6v  ' 


(II. 5-8) 


Equations  (II.4-?9a)  and  (II.4-29b),  along  with  the  above,  lead  to 


S'fp  p 

s  Ip'  q' 

s 

s' 

t  ' 
t'. 

(a) 

6v 

-  n  f1™6  r 1/2  - 

*  nv(8kBT)  °l 

p  q 

,p‘  q' 

s 

s' 

t 

t\ 

1 

(a) 

6v 

(II. 5-9) 

and 

c'i(p  q 

s 

t  ! 

(a) 

/^Svv-1/2  It 

'  nv  8kgT^  0 

p  q 

s 

t  ’ 

(a) 

s  Ip'  q’ 

s' 

t'J 

<5v 

p'  q' 

s 1 

t’. 

6 ' j 

(II. 5-10) 


Equations  (II. 5-5)  and  (II. 5-6)  reflect  the  fact  that  they 

have  been  derived  for  a  single  component  gas  in  that  they 

(k) 

involve  the  quantities,  I' ,  ,  ,  ,  ,  characteristic  of 

Jtn x,  n  )  pbp  b 

such  a  system.  The  generalization  of  (II. 5-5)  and  (1 1. 5-6) 

15 


to  a  binary  gas  mixture  requires  use  of  the  Talmi  transformation 
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as  presented  in  Ref.  12.  The  resultant  expressions  are 


[p  q  s  t 
Ip'  q'  s'  t' 


-  (-l)‘,ta+p'y  (-l)ka2(k);!(kq'q)>/2{j,  J'  ak} 

*  I  *£n;’n‘  ;psp'  s'  ^a6’av^°'  [t*  2-  5')k“ 


(II. 5-11) 


,i  f P  q  s  t  (a) 
Ip'  q'  S'  t'Jdv 


(.|)1ta4P'j  (-l)ki.2(k)si(kq'q)>/2  K  k)  l  Sl(kU'),/2 
k  [P  p  in 

I '  n ' 

(-l)e‘  (ajl)4n‘+2£'-2s'-p  q'  rf  th’ 


ini'n' ;psp  s 


(II. 5-12) 


"V  i1'2  ,(k) 


with  “«  =  (2k‘BT)  •  and  Ilni*n1;psp*s'(V“v)  as  defined 

in  Ref.  12. 

1  6v  „  t  1  6v 

In  the  above,  o' [£l  [,]  k  and  o"[£.  J.  J.Jk 

are  generalizations  of  the  collision  integrals  given  in  equations 
(1.4-15)  and  (1.6-9)  and  are  similar  except  that  in  these 
quantities,  collision  partner,  a,  is  of  species  6  ,  while  b  is 
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of  species  v  .  Thus,  the  generalization  of  the  single  component 
development  of  Part  I  to  a  binary  mixture  is  accomplished  by 
properly  considering  the  transformation  to  center  of  mass  and 
relative  coordinates  using  the  Talmi  transformation. 

Final ly. 


S' 


P  q  s  t  '(a) 
p'  q’  s'  t'Jiv 


nv(Sr?)1/5(-,)<IMtp,E  <-<)k„2(k)a(kqq')'/2K  J‘  ^ 
B  k 


zn 

z'n' 


Z  q  n  t 
z‘  q'  n*  t'jk 


(II. 5-13) 


and 


S" 


P  q  s  t  Ha) 
Ip*  q'  s'  t'J  6v 


\'Sf)1/2|-’!q+a‘P'l  (-Dk»2(k)q(kqq')'/kK  J'  k) 


X  l  £l(k?.Z')'/2(-l); 

zn 


(%4n'+2z'-2s'-p'(k)  (a  .«  Jo-ljj,  ", \ j *V  • 

‘znz'n'jpsp's’  5  v  Uq  "  tJk 

(II. 5-14) 


The  above  quantities  may,  in  principle,  be  computed,  and 
consequently,  so  also  the  transport  properties  of  a  binary  gas 
mixture  in  an  applied  field. 
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1 1. 6  Summary 


In  summary,  the  linearized  Waldmann-Snider  equation  is 
generalized  in  section  II. 1  to  a  binary  gas  mixture  in  an 
applied  magnetic  field.  The  detailed  nature  of  this  general ization 
is  examined  in  section  11.2,  and  tensor  equations  are  obtained 
for  the  perturbation  operators,  ,  appearing  in  the  perturbation 
expansion  of  the  singlet  distribution  function-density  operators, 

fk 

obtain  tensor  equations  for  the  transport  properties  of  a  binary 
gas  mixture  in  an  applied  field.  In  section  II. 4,  the  tensor 
equations  for  the  transport  properties  are  transformed  into 
scalar  equations  involving  the  quantities  S'  jP, 
and  S" 


i  ^  .  The  expressions  for  are  used  in  section  1 1.3  to 


)  lpl  q'  s'  t'J< 

p  q  s  t  (a)  .  /\n(j  finally,  in  section  II. 5,  the 

J  0  V 


Ip'  q'  S'  t' 

relative  momentum  collision  integrals  of  Part  I  are  generalized 
to  collisions  involving  two  molecules  of  possibly  different 
species.  It  is  then  possible  to  express  the  scalars, 

S! 


p  q  s  t  1(a 
IP'  q'  s'  t 1  j tS v 

the  binary  mixture  collision  integrals,  o’ 


and  S"(p’  q'  s'  •  in  U‘rmS  °f 


i  q  n  t 
f.'  q'  n'  t'J 


6v 

k 


and 


„  (l  q  n  t  )6v 
0  q*  n*  t'Jk  ' 

With  the  results  of  sections  II. 4  and  II. 5,  it  is  now 
possible,  in  principle,  to  calculate  the  scalar  quantities. 


S' 


p  q  s 

tHa) 

and  S”fp, 

q  s  t 

(a) 

p'  q'  s' 

'  5v 

IP 

q'  s'  t\ 

6\) 

,  for  a  binary 


mixture  and  using  them,  solve  a  set  of  matrix  equations  to 


obtain  the  shear  viscosity,  bulk  viscosity,  multi-component 
diffusion,  thermal  diffusion,  uniform  composition  thermal 
conductivity,  and  flow  birefringence  of  a  binary  gas  mixture 
in  the  presence  of  an  applied  magnetic  field.  The  field  effects 
on  the  shear  viscosity,  thermal  conductivity,  and  flow  birefringence 
are  the  most  widely  studied  and  reported  experimentally. 

The  uniform  composition  thermal  conductivity  ran  be  related 

to  the  experimentally  measured  thermal  conductivity  through  the 

3  9 

Stefan-Maxwel 1  equation.  ’  To  calculate  the  uniform  composition 
thermal  conductivity  to  the  same  accuracy  as  the  shear  viscosity 
and  flow  birefringence  requires  a  larger  basis  set,  due  to  the 
greater  importance  of  internal  polarizations.  Consequently, 
calculation  of  the  thermal  conductivity  requires  the  solution 
of  a  larger  set  of  scalar  equations  than  does  calculation  of 
either  the  shear  viscosity  or  flow  birefringence. 

Calculations  of  comparable  «nxuracy  can  be  performed  for 
either  the  shear  viscosity  or  flow  birefringence  using  basis 
sets  of  essentially  the  same  size  Since  the  field  effects  on 
the  shear  viscosity  have  been  widely  studied  or. pci  imentally, 
the  shear  viscosity  is  treated  in  detail  in  Part  III.  The 
discussion  there  is  intended  as  an  example  of  the  type  of 
treatment  that  is  now  possible  for  any  of  the  transport 
properties  discussed  in  Part  II. 
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Appendix  II. A  Demonstration  of  the  Relation  Expressed  in 
Equation  (II. 3-19) 

The  relation  given  in  (II. 3-19)  can  be  demonstrated  by 
means  analogous  to  those  used  in  section  7.4.c  of  Ref.  9. 

A  sketch  of  the  procedure  is  given  here. 

Using  (II. 2-11) 


2knT  ,n,  H!  ,  <H!>. 

I  ~  Vr*  \ 


*  l  Tr 
kt  ' 


Ofk  fk0) 


iLA' 


■ki 


it  dBk 
(II.A-la) 


The  RHS  of  (II.A-la)  can  be  written 


f<°>!  i  (is 


.  L  HI  )A 
kr  r  kr  -r 


-kt 


(II.A-lb) 


where  fikrAr  =  ^rAk  +  6\r^r  •  In  the  above. 


Trk  f  %  fk0)[fikr»rl?k«  '  de 

*  -(2,r'('2  Trk  f  dpk  f(0,(T.-r  |  dSr  f<°>[{ 


dp: 


x  <"kr9krlt!"kr3kr>(?kt?r><“kr3kr|6(E)tt|“kr3kr>s(Sk+Br-6ii-5;» 

*  2iT(<‘'krSkrltl“kr3kr>(*k*»r)  '  (V^VSkr^KrSk^tt^l- 


(II.A-2) 
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where  the  primes  on  ano  indicate  a  functional  dependence 

on  and  p'r  .  respectively. 

Switching  indices  in  the  term  quadratic  in  t  and 
" symmetrizing" , 


^  |  dEk  •  iz 

■  'Vr  |f  dekdprf  j"  fj 

*  <vtr3krltt|kkr3k>(A^ArKvkrgltr|ts(£)|.krg|;>«(E|t*Er-Ej;-pp 


c;  -c’ 

-kt  -r*. 


(0M0)  <* 


f'w'f 

a.  k  r 
+ 


c.  +c  . 

-Ki'  -rt 

2. 


<llkr?kr' t'  ^kr9kr^-k+-r'J 


C,  ]  \ 

+  tfV-r^kHkr'^^kr-kr'*  "~~2 


J  ^ 


-  -U*)‘h-  TrrTrk  if  %<ir  i  ifj  <lDkf(01ffJ0r<Ak*Sf) 


*  ^kArl'SrSj^  ^KrSkr' 5<E >‘h »kAr> 

j  f(0)^(0)  Q  +£ 

»  A(ek+Er-ek-s;>  *-W- 


c,  „+c  . 

-■■■■•>— —  < ,  g.  t  '  c,  >  ;A.+A  ) 
c  kr^kr  kr*kr  -k  ~r' 


(II.A-3) 


that 


where  use  has  been  made  of  the  property^ 

-  -  <-kr3ir I 1 Mkr9k  >  •  W.A-4) 

By  removing  the  "symmetrization"  from  A^+Ar  .  and  utilizing 
energy  conservation, 

Trk  l  dBk<°)C<Wrl?k«  ■  h  “  TV  f  dp,f‘0,[»vk9k  A  •  d„ 

(II.A-5) 

It  can  also  be  shown  by  similar  but  more  straightforward  means 

that 

Trk  f  "ekCtv^ At  •  dk  "  Trr  J  der40)[ VkWV  dt . 

(II.A-6) 

Equations  (II.A-5)  and  (II.A-6)  allow  the  demonstration  to 


be  completed.  Beginning  with  the  RHS  of  equation  (II.A-1), 
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PART  III. 

SHEAR  VISCOSITY  IN  A  MAGNETIC  FIELD 


III.l  Solution  of  the  Transport  Equations  for  a  Binary  Mixture- 
Example:  The  Shear  Viscosity 

The  equations  describing  the  effect  of  a  magnetic  field 
on  the  transport  properties  of  a  binary  gaseous  mixture, 
developed  in  Part  II,  are  valid  for  mixtures  in  which  both 
components  are  diatomic.  Since  a  calculation  for  a  diatom- 
diatom  mixture  (or  pure  diatomic  gas)  is  considerably  more 
difficult  than  an  atom-diatom  mixture  in  which  the  diatomic 
species  is  present  in  low  concentration,  the  solution  of 
equations  (II. 4-31)  and  (II. 4-41)  for  the  shear  viscosity 
tensor  in  this  special  case  is  considered  in  further  detail. 
Thus,  the  equations  of  the  previous  chapter  are  reduced  by 
considering  the  limit  in  which  one  of  the  diatomic  species 
becomes  spherical  while  the  other  diatomic  is  present  in 


low  concentration. 
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In  section  III. 2  an  expression  for  the  shear  viscosity  is 
obtained  subject  to  the  above  restrictions.  The  truncation  of 
the  basis  set  is  that  used  by  Hunter  1  in  a  discussion  of  a 
single  component  diatomic  gas.  In  section  III. 3,  the  collision 
integrals  involved  in  the  shear  viscosity  calculation  are 
examined  in  some  detail.  Their  calculation  in  the  classical 
limit  is  outlined  in  section  III. 4.  Finally,  a  qualitative 
comparison  of  the  field  effects  predicted  by  the  expression 
obtained  in  section  III. 2  with  experimental  observations  is 
discussed  in  section  I II. 5. 


*'  V 
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III. 2  Spherical  Components  of  the  Shear  Viscosity  Tensor 

In  this  section,  expressions  for  the  components  of  the 
shear  viscosity  of  a  binary  mixture  in  an  applied  magnetic  field 
are  obtained  using  a  truncated  basis  set.  The  results  are 
expressed  as  ratios  of  determinants.  For  the  special  case  in 
which  only  one  component  of  the  gas  is  diatomic,  while  the  other 
is  atomic,  simplifications  are  possible,  expressions  for  this 
special  case  are  given  in  detail. 

Recalling  equation  (II. 4-41),  r,  is  determined  by  the 
expansion  coefficients,  .  and  B(2^m-2  ’  where  the 

subscripts,  1  and  2,  label  the  two  species  present  in  the 
binary  mixture.  If  the  basis  set  for  each  species  is  truncated 
to  include  only  the  terms  pqst  =  2000  and  0200  ,  (requiring 
a  =  2)  equation  (II. 4-31)  yields  the  following  set  of  equations: 


♦  B(2V;i!/''??S§>!ftS"(S?8S)!f*s,(258§)!jr1  *  b^;2s"(§oSS,(2> 

(III. 2-1 ) 

J5  =  Bt°00  ciw2300,{2)  q  2  0  0  0  j  c,/  2300,(2)  2000,(2)  .2000,(2)1 

B(2)m;  r  I2000J21  B(2)m;2j_  ^2000 )22  S  (200o)22+s  (2ooo)21J 

+  B0200  c../ 2000,(2)  0200  I'  ,.2000.(2)  2000  (2)  2000  (2r 

B(2)m;r  ^0200)21  +  °(2)m;2!_  (o?.oo)22+S  (o?oc)22+s  (0200)3]  . 


(III. 2-2) 


138 


0  =  8 


2000 
( 2)m; 1 


+  B 


0200 

( 2  )m ;  1 


■  ,0200.(2)  0200. (2)  .0200.(2)]  2000  ,0200.(2) 

1  '2000/]']  +->  V2000/]]  '2000/]2  +  “(2)m;2^  '2000/]2 

S'(0200)]]^+S"{0200)j]^+S'(0200  )]2^  +  J 


,  „0200  -,,,0200.(2) 

+  B(2)m;2S  (o20o)12  • 


(III. 2-3) 


and 


0  -  ■c5;il1s‘(S5S84V*B!l)-52[s,tS588)‘|,*s“(ISgg)t|>+s-(iigJ)<f>' 


+  B, 


0200  0200  (2)  0200  "  , 0200 . (2)  , . 0200  ( 2)  , 0200  (2) 
(2)m;lb  1 0200 i21  +  B(2)m;2-  l°?00'22  S  (°200'22  S  (o20o)21 


+  im 


im  (ojl) 2 


(III. 2-4) 


r2000  and  B?°?°  can  then  be  expressed  as  ratios  of 

D(2)m;l  (2)m;2 

determinants: 


2000 

B(2)m;l  = 


ft 

A("0 

"12 

A(m) 

*13 

A(m) 

"14 

ft 

«(">) 

"22 

A(m) 

"23 

A(m) 

h24 

0 

A(m) 

h32 

A(m) 

"33 

M 

"34 

0 

A(m) 

"42 

A(m) 

rt43 

A(m) 

"44 

i*S,i 


(III. 2-5) 


in  which 


i~  2000,(2}  ,2000,(2}  2000  (2 V 

IS  (2000  )}^  +  S"(200o)5i  +S'  (2000)},' 


11 


'12 


,(III.2-7a) 


cu/2000  (2) 

v 2  0 0  0  '12 


(111. 2-7b) 


A 


(m) 

13 


13  0  0 


,  . - ,(2)  2000, (2}  2000  (2} 

S  (0: og )j^+S"{ 0200) j^'+S'  (0200) 


,(III.2-7c) 


A 

A 

A 


(m) 

14 


(m) 

21 


(m) 

22 


200P 

S"(  0200^2  ’ 


s"(I?SS)^>  . 


(III. 2-7d ) 

( III .2-7e) 

(III. 2-7  f ) 
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A("0 

"23 

= 

s"(8$83)<?>  . 

(Ill .2-7g) 

A(m) 

"24 

= 

s-(g?S8)«l>+s-(l2SS)(|>+s-(iSSg)(f)]  , 

(III. 2-7  h) 

A(m) 

"31 

= 

*  — ? 

_s-<Si8S)<f>«-(SSS8)1tfUs-(SiSg)f|>J. 

(III. 2-71) 

A(m) 

ft32 

= 

S"(:5oo){!)  . 

(III- 2-7 j ) 

A(m) 

"33 

= 

1, 0200.(2)  .0200.(2)  ,.,.0200,(2)  ■  , 
IS  (0200)|-j  +S  ( 0200 )  ^  -j  +S  (0200)|2  +  lm( 

1 

"lV 

> 

( l I I .2-7k) 

A(m) 

"34 

= 

s-(8188)<!>  , 

(III. 2-70 

A(m) 

M41 

= 

S"(2000)^  , 

(III. 2-7m) 

A(m) 

"42 

= 

[“s- (ll88)il)+s»dS88)'|,-s- (?S88)‘f ’]  . 

(III .2-7n) 

A(m) 

"43 

= 

s"<8i88)‘?>  . 

(III. 2-7o) 

A(m) 

«44 

= 

“l )  2J 

• 

(III. 2-7p) 
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Following  Ref.  2,  can  be  written  as  the  ratio  of 
determinants: 


H(m) 

H11 

M(m) 

1 2 

u(m) 

”13 

u(m) 

”14 

Hi 

n 

u(m) 

H21 

H(m) 

22 

M 

”23 

H(m) 

”24 

Hi 

n 

H(m) 

”31 

H(m) 

H32 

H(m) 

H33 

u(m) 

”34 

0 

(III. 2-8) 

w(m) 

H41 

H<m) 

”42 

H(nO 

”43 

H(m) 

H44 

0 

m 

n 

Hz 

n 

0 

0 

0 

in  which 


If  one  of  the  constituents  of  the  binary  mixture  is  atomic, 
the  term  with  pqst  =  0200  is  no  longer  an  allowed  member  of  its 
basis  set,  and  (III. 2-7)  no  longer  applies.  For  the  case  in 
which  species  1  is  atomic: 
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M(m) 
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H(m) 

"14 

*1 

M 
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u(«n) 
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u(m) 
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x2 

H(m) 
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H(m) 

"42 

44 

0 
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"1 2 

H(m) 

"14 

„(m) 

"21 
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H(m) 

"24 

— 

M(m) 

"41 

H(m) 

"42 

where  x.  = 


n. 

T 

1  "  n 


is  the  mole  fraction  of  species  i 


Letting 

c(P  q  s  t  x (a ) 
Mp*  q'  s'  t,;i 


(III. 2-10) 


r.-fP  q  s  t.\^  +  c»fP  q  s  t  \ (a n 
|_S  {p’  q’  s'  t* Mi  *  S  V  q*  S'  t  Mi  J 
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> 

Cin.2-11) 


.  n  ,^Lrl/2  0(P  q  5  r  )ia; 
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and  expanding  (I I 1.2-10) 


\  m\ 
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Now,  for  x:>  small,  equation  (III. 2-12)  can  be  expanded 
in  powers  of  x2  ,  yielding 
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Evaluation  of  A/t  is  immediate  by  inspection: 
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Defining  two  new  quantities, 
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( III .2-21c) 


For  concentrations  of  the  diatonic  species  sufficiently 
small,  only  terms  to  first  order  in  x2  need  be  kept,  i.e.,  only 
the  first  two  terms  in  (III. 2-19).  Under  these  conditions. 


-21a) 


4 


nm  depends  only  on  cross  sections  for  atom-atom  and  atom-diatom 
collisions.  The  primary  contribution  to  nm  is,  of  course, 

.  The  effects  of  the  field  and 

the  corrections  due  to  the  presence  of  the  diatomic  occur  in 
the  term  first  order  in  x2  ,  with  all  of  the  field  dependence 
isolated  in  the  iimj>  and  im$  terns. 
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III. 3  Expressions  for  the  Scalars  Involved  in  a  Computation 
of  the  Shear  Viscosity  Tensor  of  an  Atom-Diatom  Mixture 


To  calculate  the  shear  viscosity  tensor  for  a  small 
concentration  of  diatomic  species  in  a  predominantly  atomic 
gas,  the  first  two  terms  of  equation  (III. 2-19)  are  sufficient. 
These  terms  involve  ten  scalar  quantities,  one  for  atom-atom 
collisions  and  nine  for  atom-diatom  collisions.  Of  the  nine 
atom-diatom  scalars,  five  are  of  the  o’  type,  while  four  are 

of  the  a"  variety.  The  scalar  quantities  needed  for  such  a 

,  ,  .  ,2000.(2)  ,,2000.(2)  , ,2000.(2) 
calculation  are  o{2000)|  '  ,  o'(2000)^/  ,  a'(20oo)j2'  , 

,.2000.(2)  .  .0200.(2)  ,  .0200,(2)  2000  (2) 
a'(0200)^;  ,  o'(2000)^;  ,  o'(0200)^;  ,  o"(2000)^  , 

,.,2000  /O'  2000  (2)  ,0200  (2) 

o' (2000) j:/  ,  a  (0200)^  .  c’  (2000)^;  ,  where  once  again, 

1  is  the  atomic  species  and  2  the  diatomic. 

Equations  (II. 5-11)  and  (II. 5-12)  give  expressions  for 
the  above  quantities  in  terms  of  the  collision  integrals, 

0  't*  q*  n’  t’ 'k  nd  0  U*  q’  n’  t’'k  *  These  exPress10ns 
involve  the  quantities,  ^nt'n’ ipsp’s’^e^v^  *  resu1tin9 
from  the  Talmi  transformation,  six  of  which  are  needed  for 
the  calculation  discussed  here: 
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6tn,106t'rf,10 
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Using  the  above  and  the  relation 


q  n  t.6v  _  q  n  t.fiv  .  , ,1  q  n  tv 

£ '  o  n'  o'k  “  0  '£'  o  n'  o*k  6kq°  o  n'  o'k 


which  requires  k  be  zero  if  q=q’=0,  equations  (II. 5-11)  and 
(II. 5-12)  can  be  evaluated  for  the  scalars  of  interest: 
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and 
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0200.(2) 

2000^21 


m? 
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, ,0200.21 

0  ' 2  00  0 '2 


(III  .3-2.) ) 


In  (III.3-2a) ,  (III.3-2b),  (III. 3-2c ) ,  (III. 3-2g ) ,  and  (111.3-Zh), 
use  has  been  made  of  the  fact  that 


, . uuuu. < 

0  (oooo). 


=  o"(Q000°0)0y  =  0 


Equations  (1.4-15)  and  (1.6-9)  can  now  be  used  to  obtain 
expressions  for  the  nine  collision  integrals  appearing  in 
(III. 3-2)  above. 

In  general,  i.e.,  if  both  collision  partners  are  diatomic, 
equation  (1.4-15)  yields 


o'dSooC  =  TTirirr  l  -»(**’){« ... -w2  l  K^'Xp,.,.) 

V  B1  xx 1  Ja,jb'  JaJ  b 


X  x'  2 

0  0  oJ 
(-lfh’lM3 


f  ^  \\  dVSb  s*(VVx'lsasb> 
l\V 


h  M3  -h  M3  OJ  j  [y2  -  kgT  K]  S(ja'  jb'*ISdSb)]l6v  ■ 

(III. 3-4) 


k, 
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If  one  of  the  collision  partners  is  spherical,  the  reduced 
scattering  matrix,  S  ,  is  independent  of  the  angles  associated 
with  the  orientation  of  that  collision  partner,  and  the  j 
quantum  number  associated  with  rotational  states  of  that 
collision  partner  vanishes.  Now,  utilizing  the  fact  that 

^W^aV  =  exp  E2iH<W'->SaSb)]  * 

(III. 3-5) 


where  H( ja  *  jb ’ A I saSb^  is  the  9enera1i?ed  phase  shift 
introduced  by  C.  F.  Curtiss,  ^  the  o’^ooo)^  collision 
integrals  appearing  in  (III. 3-2)  can  be  written. 


2000  91  j  |,  2  _  , 

o*  (2000  )q*  =  {t—Ij  [iu(XX1){«aaI  -  (8n 2)-1  l  [o(XX  1  )  (p  ,  ) 

[>  ^  1  3 


(0  o’  0;  j  d^3  e‘^  j  dSa  exp[-2iH*1(ja(x',Sa)] 
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-hMa)\  x 


2) 


*  [y2  .  n;FTK(a)J  exp  C2iH21(ja-x|Sa)]]}  ,  ( 1 1 1 .3-6a) 

'■  B 


The  cr'(iooo)gV  collision  integrals  are  treated  in  a 
similar  manner  to  yield 


o'fiooojp1  =  I  l  a(XX' )  ji  6  -  (8it2)-1  l  [ot(XX')  (p  - , ) 

u  J  U21kb'  xx*  aa  j.  Ja 
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;  jff ’M<a)  -h"1.'4a)  o)| 

^2-2nrrK(a))1/2]  exPt2iH?.i(j;xisa)] }  (m.3-7a) 
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and 
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j  dYY2  e'Y~  j  dSb  exp[-2iH*2( j^A 1 | S^)] 
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[{Y2-2i7FTK(b)),/2'i  exp[2iH12(j£>x‘|Sb)3  }  .  (III. 3-7 b) 


The  remaining  collision  integrals  appearing  in  (III. 3-2)  are 
treated  individually  and  yield  somewhat  more  complicated 
expressions  than  those  of  (III. 3-6)  and  (III. 3-7). 
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Using  equations  (1.4-15)  and  (1.6-9),  o' (0200  and 
.2000  12 

<*"(0200)2  can  be  written 
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0200  0200 

while  a' (2000)0  and  o'(o?oo)q  become 
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From  equations  (III. 3-6a )  and  (III. 3-6b)  and  equations 
(III. 3-7a )  and  (III.3-7b)  it  is  apparent  that 
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2000  21 
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2000  12 
(  2000  )  g 


(III. 3-9) 


and 
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1000  21 
(lOOO)g 
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(III. 3-10) 


Thus,  seven  collision  integrals  need  he  evaluated  to  obtain  the 
five  spherical  components  of  the  viscosity  tensor  for  an 


atom-diatom  mixture  subject  to  an  applied  field.  In  the  next 

2000  11  2000  21 

section,  these  colMsion  integrals,  o'(2ooo)g  »  o'(20oo)g  . 

1000  21  2000  21  ..,2000,12  .,0200,21 

o  (iooo)q  .  o  (0200)2  ’  0  '0200'2  ’  0  '2000'2  ’  and 

0200  21 

o' ( 0200 )g  »  are  written  in  forms  appropriate  for  computational 
purposes. 
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III. 4  Classical  Limit  Expressions  for  the  Collision  Integrals 

Expressions  for  the  classical  limits  of  the  seven  collision 

integrals  listed  at  the  end  of  section  1 1 1. 3  are  next  considered. 

The  procedure  for  obtaining  such  classical  limits  is  discussed 

5  4 

in  detail  in  the  work  of  K.  Squire  and  R.  Wood,  and 
consequently,  is  only  sketched  here. 

An  expression  of  the  form 


[  dS  dS  e~2iH*( jajjjA  A  e2 i H ( ja j |jX |  SaSb)  , 

J  a  b 

where  A  is  any  operator  in  the  space  of  the  six  Euler  angles 
indicated  by  Sa  and  $b  ,  may  be  rewritten  as 


ff  <V5b 


e-.?iH*'jajbx|SaSb)e2i[H*(jajbA|SaSb)-H*(jajb^|SaSb)] 


X  A  e2iH(jajb>.|SaSb) 


(III. 4-1) 


It  is  shown  in  Ref.  5  that  there  exists  a  unitary  operator, 
P(A;S  S.)  ,  such  that 

a  D 


2*H{JaVlSaV 


•  PU-.Vb> e 


2ii 


(III. 4-2) 


where  n^  is  the  spherical  phase  shift.  This  makes  it  possible 
to  write  (III. 4-1)  as 
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(III. 4-3) 


in  which  A  =  P~  l(A;S  S.  )  A  P(*;S  S  );  f.  =  a'-\  ;  and  the 

a  D  3D 

generalized  angle  of  deflection,  x-(j  Jk'IS.S.)  ,  is  defined 

s  3D  3D 

such  that 


<  P(x;SaSb) 


(III. 4-4) 


and  consequently. 


- 1  2i [H*( j  j  A | S aS .)-«*( j _jbA+C | S  S  )] 

P  \x;S&Sb)e  a  b  a  b  a  b  a  b 


-iCx-(jajbx|SaSb) 


(III. 4-5) 


Next,  the  operator,  A  ,  is  replaced  by  the  classical  limit, 
A(S_S. )  ,  which  is  a  function,  ^  and  xc  ’s  replaced  by  its 
limit,  which  is  a  function  independent  of  £  .  Finally,  infinite 
sums  over  discrete  indices  are  replaced  by  corresponding  integrals 
over  continuous  variables.  Details  of  the  above  procedure  are 
discussed  in  Ref.  5. 
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Carrying  out  the  above,  the  classical  limits  of 

,  .2  000.11  .,2000,21  .  ,.  1000,21 

0  ( 2 o o o ) o  .  o  (2000)0  *  and  0  (iooo)q  are  given  by  the 

following  expressions: 
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(III  .4-6b) 


|  dye"Y  f  bdb  jY5-2YM  |  dtrote 
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(8^W  f  dS  |y2-  .-\—  K(a)1}/" 
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n  4d):  lll/?  kh 
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(III -4-6c ) 


in  which  b  ii  the  classical  impact  parameter,  ^  f  ^  is  the 

5 

dimensionless  rotational  energy  of  the  diatomic  molecule, 

X]  is  the  spherical  angle  of  deflection  for  colliding  atoms, 
X2l(Sd)  is  the  generalized  angle  of  deflection  for  an  atom- 
diatom  col  1  ision,  and  and  are  functions  defined 

in  Ref.  5. 
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2000  21  £300  12 

The  classical  limits  of  a*(o20o)2  »  ^"(0200)2  * 

0200  21  0200-21 

o' (2000)2  ,  and  a'(0200)g  are  somewhat  more  difficult  to  obta 

They  are  found  by  replacing  by  by  , 

K(i)  by  K(a)  ^  j  (2, -a)  by  lia)(2,-a),  and  D2  (Sj  by 

a  1  p  a 

D^(Sa)  »  and  again  introducing  the  generalized  anqle  of 
deflection.  The  quantities  ,  and  are  those 

of  References  4  and  5,  while  1^(2, -a)  and  D2  (S  )  are 

Ct  p  a 

the  classical  limits  of  P  _1I  (2,-a)P  and  P“1D/  'S  )P  , 

a  nK  J 

respectively.  The  classical  limits  of  the  Rg^  Wang  Chang- 
Uhlenbeck  polynomials  are  independent  of  the  arguments.  They 

are  <— hT-)Z-?A-r"')3  •  wM,e  :[J]U))ja 

is  rj,(j+l)  [4j-{ j.+i )-3]  .  Following  Reference  5  and 
o  a  a  a  a 

2000  21 

utilizing  the  above,  the  classical  limits  of  o‘(o20o)2  . 

2000- 12  ,.0200.21  .  ,.3200.21  .  .  .  .  , 

o  (o20o)2  .  0(2000)2  1  and  o'(o20o)^  may  be  ootained. 

Once  the  classical  limits  of  the  seven  collision  integrals 

discussed  in  this  section  have  been  obtained,  it  becomes 


in. 


possible  to  compute  them  numerically.  The  methods  used  by 
R.  Wood  for  a  Lennard-Jones  6-12  potential  with  P2  anisotropic 
attractive  and  repulsive  terms  may  be  extended  in  a  straightforward 
manner  to  the  collision  integrals  involved  in  a  calculation  of 
the  shear  viscosity  tensor  of  an  atom-diatom  system  subjected 
to  an  applied  magnetic  field.  Though  straightforward,  the 
computations  indicated  here  are  expected  to  be  quite  lengthy. 


1 65 


It  is  possible,  however,  to  discuss  the  nature  of  the  effect 
of  an  applied  field  on  the  viscosity  tensor  obtained  in 
section  II 1. 2  without,  carrying  out  such  a  calculation. 

This  is  discussed  in  the  following  section. 
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III. 5  Relationship  of  nm  to  Previously  Calculated  and 
Experimentally  Measured  Results 


In  this  section,  attention  is  focussed  on  the  qualitative 
behavior  of  nm  as  given  by  equation  (III. 2-19).  Both  low  and 
high  magnetic  field  limits  are  examined,  as  well  as  the  spherical 
limit  and  the  experimentally  measureable  quantities.  Qualitatively, 

4 

nm  is  shown  to  be  consistent  with  recent  computations  of  n 
in  the  absence  of  a  field,  and  with  the  measurements  of  Beenakker 
and  others  7  of  viscosities  of  gases  with  diamagnetic  diatomic 
components  in  the  presence  of  an  applied  magnetic  field. 

Recalling  equations  (III. 2-19)  and  (III. 2-20), 


=*  _atom  ra+bO  +  imiM-n  /tII  c  1X 

nm  -  n  ti  ♦  x;l-vnv^.aJJ  .  (ill. 5-1) 


■KUtkrJ  V,  , 
atom  _  ,  l  8  x  2  1 

n  -  (—« — )  — — - 


(III. 5-2) 


and  in  which  only  terms  to  first  order  in  the  mole  fraction  of 
the  diatomic  species  have  been  kept. 

In  the  limit  that  the  magnetic  field,  ,  goes  to  zero, 
tp  and  <t>  go  to  zero,  and 


(III. 5-3) 


with 


167 


no  = 


atom 

n 


{ 1  +  X  2 


(III. 5-4) 


If  expanded  in  orders  of  the  anisotropy  of  the  atom-diatom 
potential,  the  zeroth  and  first  order  terms  of  no  are, 
respectively. 
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and 


no  ( ^ )  =  0  • 


(III. 5-6) 


In  the  above,  the  collision  integrals  subscripted  with  zeros  are 
the  contribution  to  the  collision  integral  from  the  spherical 
portion  of  the  interaction  potential,  i.e.,  the  zeroth  order 
in  the  anisotropy  contribution.  Through  first  order  in  the 


r 
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anisotropy,  no  is  identical  to  n  computed  by  R.  Wood. 

(See  Appendix  III. A.}  The  second  and  higher  order  contributions 
are  expected  to  differ  slightly  due  to  the  somewhat  different 
truncation  procedure.  (See  section  1 1 1. 2.) 

In  the  spherical  limit,  i.e.,  the  limit  in  which  only 
terms  to  zeroth  order  in  the  anisotropy  of  the  potential  are 
kept,  A  .  .  .  .  ,  ,  a  •*  0  ,  and  n  1 0)  becomes 

independent  of  both  the  field  and  the  index  m  .  This  is  as 
expected,  since  the  field  effects  are  due  to  the  non-sphericity 
of  the  interaction  potential. 

In  the  large  field  limit,  nm  again  becomes  independent 
of  the  field.  For  m  =  0  ,  of  course,  it  is  identical  to  the 
low  field  limit.  However,  for  m  /  0 


hm- 


m 


natom  (1  +  x2£) 


(III. 5-7) 


Since  all  field  dependence  is  isolated  in  \ p  and  4>  , 

(III.  5-7)  can  be  written  (4/  and  $  are  both  proportional  to 

K)2 

— - —  ,  which  is  proportional  to  H^/p) 


Vo(hhw) 


atomr. 
n  il 


+  x 


/uir; 


2a(22g0oS){2)a'(i0o^)^a-(^§)(2) 


} 


(III. 5-8) 
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r 


which  is  clearly  independent  of  field.  is  ^ent^ca^ 

to  no(0)  if  the  zero  subscripts  are  removed  from  the  atom-diatom 

collision  integrals  occurring  in  no(0)  •  The  fact  that 

tWo(Hm*°0)  ^ndePendent  °T  the  field  is  explained  physically 

by  the  fact  that  the  precession  frequency  becomes  much  greater 

than  the  collision  frequency,  resulting  in  a  statistical 
8 

averaging  of  the  diatomic  species  for  large  field.  The 

net  result  is  a  saturation  of  the  Senftleben-Beenakker  effect 

7  8 

and  is  observed  experimentally.  ’ 

Turning  to  experimental  results,  the  quantities  that  have 

7  8  i  i 

been  measured  ’  are  the  changes  in  . 

and  ^(n2+n_2+fii+n_i )  induced  by  an  applied  magnetic  field. 

The  first  two  quantities  are 


1  #  \  _  4  t 

tv rir~n_2 )  ~  o'  n 


4  atom  (a+bH-bij' 


c[l+/V 


(III. 5-9a ) 


1  /  \  C 

3V n] 1 ;  -  2  n 


2  atom  (a+bi^-bi1' 


(III. 5-9b) 


c[l+4‘ 


The  interesting  point  to  note  here  is  that  since  both  4  and  4> 
are  directly  proportional  to  H^/ p  ,  the  above  quantities  are 
odd  with  respect  to  the  field.  This  is  observed  experimentally. 
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8  9 

The  third  quantity,  measured  by  Beenakker  and  co-workers,  * 
is  (to  first  order  in  x2) 

_  An  _  f4nn-(n?+r)-2+nx+n_i)|  s  x^  j [5+84>2][(a+bH2-b<H]i 

n  ‘  1  4n0  1  2  1  c[l+5*2+4*4] 

(III. 5-10) 


This  expression  is  even  in  the  field,  as  is  observed 
experimentally.  Furthermore,  in  the  low  field  limit,  —  is 
proportional  to  the  square  of  H^/p  ,  while  in  the  high  field 
limit,  it  becomes  constant: 


An  HM"~  r  (a+b)  _ b _ , 

21  C  a(1888)1(2>a*(1888)^f)a'(8188)^f)  ' 

(III. 5-11) 

(-  — )  grows  monotonical ly  from  zero  to  the  above  saturation 
n 

limit.  All  of  these  characteristics  are  consistent  with 
experimental  observations. 

It  is  interesting  to  note  that  since  both  the  high  and  low 
field  limits  are  independent  of  m  ,  except  for  m  =  0  , 


,  an,  ^  •  \(hm  *  01  -  WhM  *  "> 

(-  n  * 


(III. 5-12) 


Thus,  the  saturation  value  for  the  relative  change  in  the 
m  =  +2,  +1  spherical  components  of  the  viscosity  tensor  is 
Identical  to  the  value  measured  for  the  mixture  in  the 
experimental  apparatus  of  Beenakker  and  co-workers. 


L. 


, , , _■ ...... 


The  important  result  of  this  section  is  that  using  the 
minimum  basis  capable  of  exhibiting  field  effects,  it  has  been 
possible  to  obtain  expressions  for  the  spherical  components  of 
the  viscosity  tensor  which  behave  qualitatively  in  a  manner 
that  is  consistent  with  experimentally  observed  behavior. 

Since,  in  principle,  computations  of  the  field  effects  are 
no  more  difficult,  than  the  recent  non-field  calculations 
of  R.  Wood,  it  is  now  possible  to  obtain  theoretical  values 
for  the  field  effects  on  the  viscosity  of  gaseous  mixtures 
that  conta in  small  concentrations  of  diatomic  molecules  in 
a  predominantly  atonic  pas.  There  is  considerable  experimental 
information  with  which  these  theoretical  results  can  be 
compared. 
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1 1 1. 6  Summary 

To  summarize,  in  section  1 1 1. 2  an  expression  (equation 
(III. 2-19))  is  obtained  for  the  spherical  components  of  the 
shear  viscosity  tensor  of  an  atom-diatom  mixture  in  the 
presence  of  an  applied  magnetic  field.  This  expression  is 
shown  in  section  III. 5  to  be  qualitatively  consistent  with 
experimental  observations  of  these  field  effects.  A  calculation 
using  equation  (III. 2-19)  requires  the  consideration  of  ten 
scalars. 


/2000,(2)  i/2000\(2)  ,/2 000\(2)  i/2000\(2)  ,,0200\(2) 
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,./0200\(2)  „/2000x(2)  „,2000x(2)  h/2000\(2)  a 

1  '0200!  »  0  ' 2 000 1 21  ’  a  '20oo7-|2  *  0  '020q/]2  ’  and 


../0200*(2) 

0  v  2000  721  > 

which,  under  closer  examination  in  section  III. 3,  prove  to 
require  the  evaluation  of  seven  collision  integrals. 
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A  computation  similar  to  that  recently  completed  by 
4 

R.  Wood  requires  expressions  for  the  classical  limits  of 
the  above  collision  integrals.  The  procedure  for  obtaining 
such  classical  limit  expressions  for  the  collision  integrals 
examined  in  section  I II. 3  is  outlined  in  section  III. 4.  It 

5 

relies  heavily  upon  the  work  of  K.  Squire,  and  the  reader 
is  referred  to  that  work  for  a  more  detailed  treatment. 


The  fact  that  equation  (III. 2-19)  is  consistent  with 
the  available  experimental  results  is  most  satisfying.  It 


can  also  be  shown  that  the  classical  limit  expression  for 
o' (2000)^  leads  to  an  expression  for  natom  that  is 


consistent  with  those  of  other  treatments.  The  RHS  of  (III.4-6a) 

is  equivalent  to  ^J^is’kT  where  is  an  omega 

®  1 1 
integral  defined  in  Reference  2  and  equivalent  to  the  [i^;(s) 


of  Chapman  and  Cowl  inn 

.  atom  r’-ikBT 
that  ’1  -  \  — p - 


Using  this  result  and  the  fact 


y-r~  • the  class1ca 


1  limit 


expression  for 


atom  , 

r.  becomes 


k  T 

atom  5  B 
nCl.  =  J2^J  ' 


(III. 6-1) 


which  is  precisely  the  expression  given  in  Ref.  2. 
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Thus,  the  dominant  contribution  to  the  shear  viscosity 
tensor  is  consistent  with  other  theoretical  developments, 
while  the  qualitative  behavior  is  consistent  with  experimental 
observations.  It  remains  to  be  seen  if  a  numerical  calculation 
will  lead  to  results  that  are  numerically  close  to  those 
obtained  by  experiment. 


_ 
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Appendix  III. A  Demonstration  of  the  Equivalence  of  no(0) 
and  nW000(0) 

WOOD 

Since  there  are  no  contributions  to  no  and  n  that 
are  first  order  in  the  anisotropy,  no  and  n  are  identical 
through  first  order  in  the  anisotropy  if  no(0)  and  n^>D(0) 
are  the  same.  Equation  (I II. 5-5)  can  be  written 


n0(0)  =  n(af°m  n+x2n!(0)]  , 


(III.A-1) 


k  T 

where  n^0”'  =  ^  ~(2  2")  ^  identical  to  no  ^  ,  and 
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Equations  (III. A-4a )  and  (III. A-4b)  can  be  added  to  yield 


m(0)  =  {m2[H,(0)+2H;?(0)3-[H3(0)Hu(0)-m1m2H2(0)2]  [ntiHi(O)]'1} 
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where  the  H^(0)  are  the  spherical  limits  of  the  quantities 
4 

defined  by  R.  Wood. 

Finally, 

nt(0)  =  niWOOO(0)  ,  (III.A-6) 

and  consequently, 

n0(0)  =  nW00D(0)  .  (III.A-7) 
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CONCLUSION 

In  Part  I,  the  collision  integrals,  developed  by  Hunter  and 

Snider^  and  essential  to  Hunter's  treatment  of  a  single  component 

2 

diatomic  gas  in  an  applied  field,  are  expressed  in  terms  of  the 

3 

reduced  scattering  matrix  of  Curtiss  and  co-workers.  The 
introduction  of  several  operators  leads  to  considerable 
simplification  of  the  expressions  for  the  collision  integrals. 

The  intent  of  this  work  is  to  move  closer  to  calculations  of  the 
transport  properties  of  dilute  diatomic  gases  in  applied  fields. 
Recognizing  the  difficulties  inherent  in  carrying  out  a 
calculation  which  requires  treatment  of  diatom-diatom  interactions 
leads  to  the  aim  of  first  discussing  a  physical  situation  which 
requires  only  the  consideration  of  atom-atom  and  atom-diatom 
interactions.  This  leads  to  the  development  in  Part  II  of 
expressions  for  the  transport  properties  of  binary  gaseous 
mixtures  in  applied  fields. 

Part  II  is  an  extension  of  Hunter's  work  on  single  component 
2 

gases  to  binary  mixtures.  As  in  Hunter's  work,  scalar  equations 
for  the  shear  viscosity  and  thermal  conductivity  tensors  are 
obtained.  In  addition,  scalar  equations  are  obtained  for  the 
bulk  viscosity,  mul ti -component  and  thermal  diffusion,  and  flow 
birefringence  tensors  of  binary  mixtures  in  applied  fields.  The 
motive  for  discussing  mixtures  is  the  desire  to  treat  systems 
composed  of  a  diatomic  species  in  low  concentration  in  a 
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predominantly  atomic  gas.  Such  systems  require  the  treatment  of 
atom-diatom  and  atom-atom  interactions  only  in  a  calculation  of 
the  transport  properties.  The  shear  viscosity  of  such  a  system 
in  an  applied  field  is  examined  in  detail  in  Part  III. 

An  expression  for  the  shear  viscosity  of  an  atom-diatom 
mixture  is  obtained  in  Part  III.  As  mentioned  in  the  summary, 

the  expression  obtained  is  consistent  with  both  experimental 

4 

results  and  the  numerical  treatment  of  such  systems  in  the 

5 

absence  of  a  field  recently  completed  by  R.  Wood.  Treatments 

of  single  component  systems  and  at  least  one  treatment  of  binary 

mixtures^  have  also  yielded  results  in  qualitative  agreement 

with  experiment.  This,  however,  is  the  first  treatment  of 

2 

mixtures  using  techniques  which  require  no  truncations  prior 
to  the  tensor  analyses. 

A  sketch  is  also  given  in  Part  III  of  a  procedure  which  can 
lead  to  the  eventual  calculation  of  the  collision  integrals 
occurring  in  the  expression  for  the  shear  viscosity.  A  large 
amount  of  work,  however,  must  yet  be  done  in  applying  the 

5 

numerical  methods  of  R.  Wood  to  a  calculation  of  the  shear 
viscosity  of  an  atom-diatom  mixture  in  an  applied  magnetic  field. 
This  leads  to  a  consideration  of  possible  future  work  related 
to  the  present  development. 

Clearly,  the  next  step  in  this  development  is  the  actual 
calculation  of  the  transport  properties  of  dilute  gases  in  applied 
fields.  The  calculation  that  seems  most  immediately  feasible  is 
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that  of  the  shear  viscosity  of  atom-diatom  mixtures  in  an  applied 
magnetic  field.  Other  transport  properties  of  aton-diatom 
mixtures  are  possible  once  the  scalar  equations  of  fart  II  have 
been  solved  in  analogy  with  the  shear  viscosity  development  of 
Part  III. 

Calculation  of  the  transport  properties  of  sinjle  component 
diatomic  gases  in  applied  fields  requires  considerably  more 
complex  computer  programs  than  those  needed  for  the  atom-diatom 
mixtures.  In  principle,  however,  these  are  also  possible. 

Part  I  is  a  step  in  that  direction. 

The  Senftleben-Beenakker  effect  is  felt  to  be  a  potentially 
sensitive  probe7  of  internolecular  potentials.  Its  use  as  such, 
however,  is  limited  at  present  by  the  relatively  small  amount 
of  theoretical  work4  in  this  area.  The  present  need  is  for 
calculations  to  go  along  with  the  considerable  experimental 
information  already  available.  This  thesis  is  an  attempt  to 
bridge  some  of  the  gap  between  experimental  evidence  and 
theoretical  understanding  of  the  Senftleben-Beenakker  effect. 
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